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The purpose of this project was to explore the range of characteristics present in the spectral energy 
distributions (SEDs) of quasars. The final results are included in three published papers: 

1 ) The Far-Infrared Spectral Energy Distributions of X-ray-selected Active Galaxies, J. K. Kuraszkiewicz, 
{\bf B. J. Wilkes}, E.J. Hooper, et al., 2003, ApJ, 590, 128 

This paper reports the ISO results on hard X-ray selected AGN which are less biased against red/obscured 
objects than other selection wavebands. We find that, as predicted, the IR continuum of these sources 
extends to redder sources than in optically/radio selected sample. This indicates that the latter samples 
miss a portion of the population which is fainter in the optical but can be easily picked up in the hard X-ray. 
The range of IR SEDs is roughly consistent with reddening of the IR continuum up to column densities of 
around 10 A 23 /cm A 3. Modeling of the full SED using dusty disk models demonstrated that varying the 
viewing angle can explain the observed SEDs, though rather large disks are required to fit the cooler, long 
wavelength emission. From the fits we can obtain estimates of the mass and inclination of the system. 

2) The Far-Infrared emission of Radio Loud and Radio Quiet Quasars; Polletta, M., Courvoisier, T., J-L., 
{\bf Wilkes}, B.J. \& Hooper, E.J., AA, 362, 75 

Abstract: 

Continuum observations at radio, millimetre, infrared and soft X-ray energies are presented for a sample of 
22 quasars, consisting of flat and steep spectrum radio loud, radio intermediate and radio quiet objects. The 
primary observational distinctions, among the different kinds of quasars in the radio and IR energy domains 
are studied using large observational datasets provided by ISOPHOT on board the Infrared Space 
Observatory, by the IRAM interferometer, by the sub-millimetre array SCUBA on JCMT, and by the 
European Southern Observatory (ESO) facilities IRAC1 on the 2.2 m telescope and SEST. The spectral 
energy distributions of all quasars from radio to IR energies are analyzed and modeled with non-thermal 
and thermal spectral components. The dominant mechanism emitting in the far/mid-IR is thermal dust 
emission in all quasars, with the exception of flat spectrum radio loud quasars for which the presence of 
thermal IR emission remains rather uncertain, since it is difficult to separate it from the bright non-thermal 
component. The dust is predominantly heated by the optical/ultraviolet radiation emitted from the external 
components of the AGN. A starburst contributes to the IR emission at different levels, but always less than 
the AGN (<= 27%). The distribution of temperatures, sizes, masses, and luminosities of the emitting dust 
are independent of the quasar type. 

3) The ISO view of Palomar-Green quasars, M. Haas, U. Klaas, S.A.H. Muller, F. Bertoldi, M. Camenzind, 
R. Chini, O. Krause, D. Lemke, K. Meisenheimer, P.J. Richards, {\bf B. Wilkes}, 2003, A\&A, 402, 87 

Mining the ISO data archive we provide the complete ISO view of PG quasars containing 64 infrared 
spectral energy distributions between 5 and 200 mu m. About half of the sample was supplemented by 
MAMBO and SCUBA (sub-)millimetre data. Since the PG quasars were selected optically, the high 
infrared detection rate of more than 80% suggests that every quasar possesses luminous to hyper-luminous 
dust emission with dust masses comparable to Seyferts and ultra-luminous IR galaxies (ULIRGs). The gas 
to-dust mass ratio (of those sources where CO measurements are available in the literature) is consistent 
with the galactic value providing further evidence for the thermal nature of the IR emission of radio quiet 
quasars. The SEDs represent templates of unprecedented detail and sensitivity. The power-law like near- to 
mid-IR SEDs (Fnu - nu alpha) are smooth up to far-infrared wavelengths, favoring dust heating by the 
central AGN, and we conclude that, in particular for our hyperluminous quasars at z=l , starbursts play only 
a minor role for powering the dust emission, even in the FIR. The IR spectral slopes alpha(l-lOum) range 
from -0.9 to -2.2 with a mean of -1.3 +/- 0.3. They neither correlate with the optical spectral slope 
alpha(0.3-l um), nor with the IR luminosity, nor with the FIR/MIR luminosity ratio, nor with inclination- 
dependent extinction effects in the picture of a dusty torus. We suggest that the diversity of the SEDs 
reflects largely the evolution of the dust distribution, and we propose a classification of the SED shapes as 
well as an evolutionary scheme in which this variety can be understood. During the evolution the 
surrounding 


dust redistributes, settling more and more into a torus/disk like configuration, while the SEDs show an 
initial FIR bump, then an increasing MIR emission and a steeper near- to mid-infrared slope, both of which 
finally also decrease. Strikingly, based on the sensitive ISO data now we do not only see the coarse IR 
differences between ULIRGs and quasars, but also the details and a possible evolution of the dust 
distribution and emission even among the optically selected PG sample. Regarding cosmic evolution, our 
hyper-luminous quasars in the "local” universe at z=l do not show the hyper-luminous (LFIR >? 1013 
Lsun) starburst activity inferred for z=4 quasars detected in several (sub-)millimetre surveys. In view of 
several caveats this difference should be established further, but it already suggests that in the early dense 
universe stronger merger events led to more powerful starbursts accompanying the quasar phenomenon, 
while at later cosmic epochs any coeval starbursts obviously do not reach that high power and are outshone 
by the AGN. 


SEE ADDITIONAL HARD COPY MATERIAL ATTACHED. 



I. 


Program Objectives 


The goal of the proposal is to perform ISO spectroscopic studies, including data analysis and model- 
ing, of star formation regions using an ensemble of archival space-based data from the Infrared 
Space Observatory’s Long Wavelength Spectrometer and Short Wavelength Spectrometer, but in- 
cluding as well some other spectroscopic databases. Four kinds of regions are considered in the stud- 
ies: (1) disks around more evolved objects; (2) young, low or high mass pre-main sequence stars in 
star formation regions; (3) star formation in external, bright IR galaxies; and (4) the galactic center. 
One prime focus of the program is the OH lines in the far infrared. 

The program had the following goals: 

1) Refine the data analysis of ISO observations to obtain deeper and better SNR results on se- 
lected sources. The ISO data itself underwent “pipeline 10” reductions in early 2001, and addi- 
tional “hands-on data reduction packages” were supplied by the ISO teams in 2001. The Fabry- 
Perot database is particularly sensitive to noise and slight calibration errors. 

2) Model the atomic and molecular line shapes, in particular the OH lines, using revised monte- 
carlo techniques developed by the SWAS team at the Center for Astrophysics; 

3) Attend scientific meetings and workshops; 

4) Do E&PO activities related to infrared astrophysics and/or spectroscopy. 


II. Program Achievements 

During this period we completed extending the SWAS (Submillimeter Wave Astronomy Satellite) 
montecarlo radiative transfer code to include: (1) OH; (2) the far IR lines of H 2 O; (3) [OI]; and we’re 
in the process of adding (4) the far infrared lines of CO. Dr. Eduardo Gonzalez-Alfonso, from the 
University of Alcala, Spain, was in residence in Cambridge for the spring and summer of 2003, and 
together we made dramatic progress in the analysis of ISO spectra. As a result we are able to model 
theoretically the emission and absorption from these species in molecular clouds whose structures 
we can adjust to ascertain the physical parameters most suitable for the various data observed. Fig- 
ure 1 shows the analysis of the spectrum of Arp220, from a paper soon to appear in ApJ. We also 
completed analysis of the spectrum NGC1068, and a paper is soon to appear in ApJ. Copies of both 
these major efforts are attached; for reprints of other papers please refer to the journal web pages. 

Also during this program we spent considerable effort helping to prepare for SIRTF observations, 
including IRAC observations of extragalactic sources, and including the proposed Early Release Ob- 
servations (EROs). We have since been actively involved in reduction and analysis of Spitzer data, 
as reflected on the publications list below. 

As a result, we have fully completed with great success all of the original goals as stated in the pro- 
posal. In fact the program was much more productive than we had originally anticipated, due to the 
great success we had in developing the modeling code, and in polishing the ISO data so that we could 
squeeze the maximum out of its spectral scans. 


III. Related Publications 

“The Infrared Lines of OH: Diagnostics of Molecular Clouds Compositions in Infrared Bright Galax- 
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ies,” Smith, H.A., Gonzalez-Alfonso, E., Fischer, J., Ashby, M., Dudley, C., and Spinoglio, L., 
in Proceedings of the Workshop of the Neutral ISM in Starburst Galaxies, Marstrand, Sweden, 
2003. 

“Mid and Far infrared Spectroscopy of Seyfert Galaxies,” Spinoglio, L., Malkan, M., Smith. H.A.. 

Fischer, J., in the ASP Conference Series, Proceedings of the Paris Conference on AGN, eds: 
Cohn, Combes and Shlosman, 2003. 

“Public Attitudes Towards Space Science,” Smith. H.A.. Space Science Reviews, 105, 1-2, 493, 2003. 

“The Far-Infrared Emission Line and Continuum Spectrum of the Seyfert Galaxy NGC 1068,” Spi- 
noglio, L., Malkan, M., Smith. H.A. . Fischer, J., submitted to Ap.J (2003). 

“Transmittance of Thick Metal Meshes of Various Shapes and Thicknesses in the Infrared Wave- 
length Region,” Sternberg, O., Shah, J., Moller, J., Grebel, H., Stewart, K., Fischer, J., Reb- 
bert, M., Smith. H.A .. and Fettig, R., TDW2003: International Workshop on Thermal Detectors 
for Space Based on Planetary, Solar and Earth Science Applications, 2003. 

“The Infrared Lines of OH: Diagnostics of Molecular Clouds Compositions in Infrared Bright Galax- 
ies,” Smith. H.A. . Gonzalez-Alfonso, E., Fischer, J., Ashby, M., Dudley, C., and Spinoglio, L., 
in Proceedings of the Workshop of the Neutral ISM in Starburst Galaxies, Marstrand, Sweden, 
2003. 

“The Effects of Dust in Infrared Luminous Galaxies: An Integrated Modeling Approach,” Satyapal, 

S., Dudley, C., Fischer, J., Luhman, M., Smith. H.A .. Astrophysics of Dust, Estes Park, Colo- 
rado, May 26 - 30, Edited by Adolf N. Witt, 2003. 

“The Far Infrared Lines of OH as Molecular Cloud Diagnostics,” Smith. H.A. . Ashby, Fischer, Gon- 
zalez, Spinoglio, Dudley, The Astrochemistry of External Galaxies, 25 th meeting of the IAU, 
Joint Discussion, Sydney, Australia, 2003. 

“Early Galactic Science Results from the Infrared Array Camera for SIRTF,” Megeath, S.T., Allen, 
L.E., Carey, S.J., Deutsch, L.K., Hora, J.L., Fazio, G.G., Forrest, W.J., Marengo, M., Melnick, 

G. J., Patten, B.M., Pipher, J.L., Reach, W.T., Smith. H.A. . Stauffer, J.R., Willner, S.P., BAAS, 
203, 2211, 2003. 

“Extragalactic Science Results with IRAC,” Willner, S.P., Arendt, R., Ashby, M., Barmby, P., Eisen- 
hardt, P., Stern, D., Fazio, G.G., Forrest, W.J., Hora, J.L., Huang, J., Huchra, J.P., Im, M., 
Pahre, M., Pipher, J.L., Reach, W.T., Stauffer, J., Smith. H.A .. Surace, J., Tollestrup, E.V., 
Wang, Z., Wilson, G., Yan, L., BAAS, 203, 2210, 2003. 

“The Effects of Dustin Infrared Luminous Galaxies: An Integrated Modeling Approach,” Satyapal, 

S., Dudley, C., Fischer, J., Luhman, M., Smith. H.A. . Astrophysics of Dust, Estes Park, Colo- 
rado, May 26 - 30, 2003. Edited by Adolf N. Witt. 

“The Far Infrared Lines of OH as Molecular Cloud Diagnostics,” Smith. H.A.. Ashby, Matt, Fischer, 
Jackie, Gonzalez, Eduardo, Spinoglio, Luigi, Dudley, Chris, The Astrochemistry of External 
Galaxies, 25th meeting of the IAU, Joint Discussion 21, 23 July 2003, Sydney, Australia. 

“First Extragalactic Science from IRAC, The Infrared Array Camera on SST,” Willner, S., ... Smith. 

H. , et al., BAAS, 2004. 

“First Galactic Images from IRAC, SST’s Infrared Array Camera,” Megeath, T., ... Smith. H.A ., et al., 
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BAAS, 2004. 


“DR21: A Major Star Formation Site Revealed by Spitzer, ” Marston, A.P., Reach, W.T., Noriega- 

Crespo, A. Rho, J., Smith, HA, Melnick, G., Fazio, G., Rieke, G., Carey, S., Rebull, L., Muze- 
rolle, J., Egami, E„ Watson, D.M., Pipher, J.L., Latter, W.B., Stapelfeldt, K., ApJS 154, 333, 
2004. 

“The Off-Nuclear Starbursts in NGC 4038/4039 (The Antennae Galaxies),” Wang, Z., Fazio, G.G., 
Ashby, M.L.N., Huang, J.-S., Pahre, M.A., Smith. H.A.. Willner, S.P., Forrest, W.J, Pipher, 
J.L., and Surace J.A., ApJS 154, 193, 2004. 

“IRAC Observations of M81,” Willner, S.P., Ashby, M.L.N., Barmby, P., Fazio, G.G., Pahre M., 

Smith. HA.. Kennicutt, Robert C. Jr., Dale, Daniel A., Draine, B.T., Regan, Michael W., Mal- 
hotra, S., Thornley, Michele D., Appleton, P.N., Frayer D., Helou, G., Storrie-Lombardi, L., 
Stolovy S., ApJS 154, 222, 2004. 

“The Anatomy of Star Formation in NGC 300,” Helou, G., Roussel, H., Appleton, P., Frayer. D., 

Stolovy, S., Storrie-Lombardi, L., Hurt, R., Lowrance, P., Makovoz, D., Masci, F., Surace, J., 
Gordon, K.D., Alonso-Herrero, A., Engelbracht, C., Misselt, K., Rieke, G., Rieke, M., Willner, 

S., Pahre, M., Ashby, M., Fazio, G.G., Smith. H.A.. ApJS 154, 253, 2004. 

“The IRAC Shallow Survey,” Eisenhardt, P.R., Stern, D., Brodwin, M., Fazio, G., Rieke, G., Rieke, 

M., Werner, M., Wright, E., Allen, L., Arendt, R., Ashby, M., Barmby, P., Forrest, W., Hora, J., 
Huang, J., Huchra, J., Pahre, M., Pipher, J., Reach, W., Smith. H. A. . Stauffer, J., Wang, Z., 
Willner, S., Brown, M.J.I., Dey, A., Jannuzi, B.T., and Tiede, G.P, ApJS 154, 48, 2004. 

“The Infrared Array Camera (IRAC) for the Spitzer Space Telescope,” Fazio, G.G., Hora, J.L., Allen, 

L. E., Ashby, M.L.N., Barmby, P., Deutsch, L.K., Huang, J.-S., Kleiner, S., Marengo, M., 
Megeath, S.T., Melnick, G.J., Pahre, M.A., Patten, B.M., Polizotti, J., Smith . H.A. . Taylor, 

R.S., Wang, Z., Willner, S.P., Hoffmann, W.F., Pipher, J.L., Forrest, W.J., McMurty, C.W., 
McCreight, C.R., McKelvey, M.E., McMurray, R.E., Koch, D.G., Moseley, S.H., Arendt, R.G., 
Mentzell, J.E., Trout-Marx, C., Losch, P., Mayman, P., Eichhorn, W., Krebs, D., Jhabvala, M., 
Gezari, D.Y., Fixen, D., Flores, J., Shakoorzadeh, K., Jungo, R., Hakun, C., Workman, L., 
Karpati, G., Kichak, R., Whitley, R., Mann, S., Tollestrup, E.V., Eisenhardt, P., Stern, D., Gor- 
jian, V., Bhattacharya, B., Carey, S., Nelson, B.O., Glaccum, W.J., Lacy, M., Lowrance, P. J., 
Laine, S., Reach, W.T., Stauffer, J.R., Surace, Wilson, G., Wright, E.L., Hoffman, A., Domingo, 
G., and Cohen, M., ApJS 154, 10, 2004. 

“DR21-IRS1: Spitzer-IRAC Four-Color Images of the Origin of the Massive Outflow and its Cluster of 
Embedded Stars,” Smith. H.A .. Allen, L.E., Fazio, G., Melnick, G., Marston, A.P., Gutermuth, 
R., Pipher, J., Watson, D., Carey, S., Noriega-Crespo, A., BAAS, 204, 6110, 2004. 

“The Far-IR Spectrum of Arp 220,” Gonzalez-Alfonso, E., Smith. H.A. . Fischer, J., Cernicharo, J., 
BAAS, 204, 4020, 2004. 

“Spitzer Space Telescope Observations of the Antennae Galaxies,” Wang, Z., Fazio, G.G., Ashby, 

M. L.N., Huang, J.-S., Pahre, M.A., Smith. H.A .. Willner, S.P., Forrest, W.J., Pipher, J.L., 
Surace, J.A., BAAS, 204, 4017, 2004. 

“Mid-Infrared Observations of the Halo of NGC 5907,” Ashby, M.L.N., Pipher, J.L., Forrest, W.J., 
Stauffer, J.R., Barmby, P., Willner, S.P., Fazio, G.G., Smith. H.A ., Arendt, R.G., Bock, J.J, 
BAAS, 204, 4004, 2004. 
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“IRAC Extragalactic GTO Program: From Nearby Galaxies to the High Redshift Universe,” Pahre, 
M.A., Ashby, M.L.N., Barmby, P., Fazio, G.G., Huang, J.-S., Smith. H.A .. Wang, Z., Willner, 
S.P., Pipher, J.L., Forrest, W.J., BAAS, 204, 3310, 2004. 

“Detecting Brown Dwarfs, Disks and Protostars with Spitzer: First Galactic Results from the IRAC 
GTO Program,” Megeath, S.T., Allen, L.E., Calvet, N., Deutsch, L.K., Fazio, G.G., Hartmann, 
L., Hora, J.L., Melnick, G.J., Patten, B.M., Sicilia-Aguilar, A., Smith. H.A .. Forrest, W.J., 
Gutermuth, R.A., Peterson, D.E., Pipher, J.L., Stauffer, J.R., BAAS, 204, 3306, 2004. 

“The Far Infrared Spectrum of Arp 220,” Gonzalez-Alfonso, E., Smith. H .. Fischer, J., and Cer- 
nicharo, J., ApJ 2004 (in press). 

Patent: Method for Fabricating Metallic Meshes Infrared Optics, 2004, Smith. H.A.. Fischer, 
J., Rebbert, M., Stewart, K., Sternberg, O., and Muller, D. Navy Case #84,958 


IV. Talks and Posters (2004-5) 


“Spitzer Views of DR21” Presented at the Cores, Disks, Jets, and Outflows Meet- 

ing, Canada, 2004 

“Massive Star Formation in DR21” Presented at the Star Formation Workshop, Cambridge, 

2004 

“Spitzer IRAC Views of Star Formation” Presented at the “ESA - DUSTY Workshop on Herschel 

and ALMA,” 2004 

Also: Ten posters at American Astronomical Society meetings (see above). (See prior year reports 
for a more lengthy list of prior meetings, etc.) 
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180 0*n> 180 


Figure 1: The Observed ISO spectrum of Arp220 (continuum subtracted), and overlayed with the 
modeled lines. Excellent agreement is found across the entire range. 
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The far-infrared spectrum of Arp 220 1 
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Universidad de Alcala de Henares, Departamento de Fisica, Campus Universitario, 
E-28871 Alcala de Henares, Madrid, Spain 

eduardo . gonzalezfiuah . es 
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ABSTRACT 

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 
shows absorption in molecular lines of OH, H 2 O, CH, NH, and NH 3 , as well 
as in the [O I] 63 pm line and emission in the [C II] 158 line. We have 
modeled the continuum and the emission/absorption of all observed features by 
means of a non-local radiative transfer code. The continuum from 25 to 1300 
pm is modeled as a warm (106 K) nuclear region that is optically thick in the 

'Visiting Astronomer, Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 
02138. 

2 CSIC, IEM, Dpto. Astrofisica Molecular e Infrarroja, Serrano 123, E- 28006 Madrid, Spain. 
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far-infrared, attenuated by an extended region (size 2") that is heated mainly 
through absorption of nuclear infrared radiation. The molecular absorption in 
the nuclear region is characterized by high excitation due to the high infrared 
radiation density. The OH column densities are high toward the nucleus (2 — 
6 x 10 17 cm -2 ) and the extended region (~ 2 x 10 17 cm' 1 2 ). The H 2 0 column 
density is also high toward the nucleus (2 - 10 x 10 17 cm' 2 ) and lower in the 
extended region. The column densities in a halo that accounts for the absorption 
by the lowest lying levels are similar to what are found in the diffuse clouds 
toward the star forming regions in the Sgr B2 molecular cloud complex near the 
Galactic Center. Most notable are the high column densities found for NH and 
NH 3 toward the nucleus, with values of ~ 1.5 x 10 16 cm' 2 and ~ 3 x 10 16 cm' 2 , 
respectively, whereas the NH 2 column density is lower than ~2x 10 15 cm' 2 . A 
combination of PDRs in the extended region and hot cores with enhanced H 2 0 
photodissociation and a possible shock contribution in the nuclei may explain 
the relative column densities of OH and H 2 0, whereas the nitrogen chemistry 
may be strongly affected by cosmic ray ionization. The [C II] 158 /im line is 
well reproduced by our models and its “deficit” relative to the CII/FIR ratio in 
normal and starburst galaxies is suggested to be mainly a consequence of the 
dominant non-PDR component of far-infrared radiation, although our models 
alone cannot rule out extinction effects in the nuclei. 

Subject headings: galaxies: abundances — galaxies: individual (Arp 220) — 
galaxies: ISM — galaxies: starburst — infrared: galaxies — radiative transfer 


1. Introduction 

With a redshift of z = 0.018, Arp 220 (IC 4553/4) is the nearest and one of the best 
studied ultraluminous infrared galaxies (ULIRGs). The tails observed in the optical, together 
with the double highly-obscured and compact (0"3) nuclei observed in the near and mid- 
infrared, as well as in the millimeter, strongly suggest that the enormous luminosity of Arp 
220, ~ 10 12 L©, is the result of galactic merging. Nevertheless, the concrete physical process 
responsible is still a matter of debate: the proposed sources are hidden active nuclei and/or 
bursts of star formation. 


1 Based on observations with the Infrared Space Observatory, an ESA project with instruments funded 

by ESA Member States (especially the principal investigator countries: France, Germany, Netherlands, and 
the United Kingdom) and with the participation of ISAS and NASA. 
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Molecular observations of Arp 220 provide unique clues to the physical and chemical 
processes ocurring in the nuclei and their surroundings. In the millimeter region, CO obser- 
vations have been carried out with increasingly high angular resolution (Radford, Solomon, 
& Downes 1991; Scoville et al. 1991; Okumura et al. 1994; Scoville, Yun, & Bryant 1997; 
Sakamoto et al. 1999). In particular they have shown that, on the one hand, CO emission 
arises from a region significantly more extended than the nuclei (~ 3 — 4"), and on the other 
hand that the CO(2-l) to CO(l-O) intensity ratio is lower than 1, thus suggesting that CO 
mainly traces low density regions (< 10 3 cm -3 ). Observations of molecules with high dipole 
moment, like CS and HCN, have revealed that the fraction of molecular gas contained in 
dense clouds (n(H 2 )> 10 4 cm -3 ) is much larger than in normal galaxies, yielding ~ 10 10 M 0 
of dense gas (Solomon, Radford, & Downes 1990; Solomon, Downes, & Radford 1992). Rad- 
ford et al. (1991b) found that HCN(l-O) and HCO + (1-0) peak strongly toward the nuclei, 
but also show low-level extended emission. More recently, Aalto et al. (2002) have detected 
emission from the high density tracers HNC and CN, and the relatively low HCN/CN and 
HCN/HNC intensity ratios were attributed to widespread PDR chemistry. 

The launch of the Infrared Space Observatory (ISO) opened a new window for the study 
of the physical and chemical properties of ultraluminous infrared galaxies. Despite the lack 
of angular and spectral resolution, the observations of Arp 220’s far-infrared spectrum from 
40 to 200 pm (Fischer et al. 1997, 1999) and of a number of individual lines in the SWS 
range (Sturm et al. 1996) provided new insights in our understanding of the ionic, atomic 
and molecular content of the galaxy. These wavelength regions are of great interest, because 
the bulk of the enormous luminosity is emitted in the far-infrared, and also because they 
contain lines of interesting molecular, ionic, and atomic species. Skinner et al. (1997) 
reported the detection of the 35 pm OH line in Arp 220. Fischer et al. (1997, 1999) found 
that the far-infrared molecular absorption lines of OH, H 2 0, CH, and NH 3 are significantly 
stronger in Arp 220 than in less luminous infrared-bright galaxies while the fine structure 
lines from ionic species are, to the contrary, extremely weak. Luhman et al. (1998, 2003) 
found that, relative to the far-infrared luminosity, the [C II] in ULIRGs is typically nearly 
an order of magnitude weaker than in lower luminosity infrared-bright galaxies. Sturm et 
al. (1996) reported the detection of two ortho-H 2 pure rotational lines, indicating that high 
masses of gas are subject to PDR conditions and/or shock activity. 

The physical and chemical processes that account for the rich molecular far-infrared 
spectrum of Arp 220 can be better understood if quantitative values of the column densities 
of the above species, as well as their excitation conditions, are estimated. The presence of 
OH and, to some extent, of H 2 0, may be indicative of PDR and/or diffuse interstellar cloud 
chemistry, and their column densities potentially give an estimate of the UV field in the 
source. On the other hand, large amounts of H 2 0 are produced in non-dissociative shocks 
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(e.g. Cernicharo et al. 1999), where the OH abundance is also enhanced (Watson et al. 1985). 
The OH abundance is expected to be generally higher than that of H 2 0 in fast dissociative 
shocks (Neufeld & Dalgarno 1989). H 2 0 ice in grain mantles may also efficiently return 
to the gas phase through sublimation of mantles in “hot core” regions. But whatever the 
characteristics of the regions producing the observed molecular features, the lines under study 
lie at wavelengths where the enormous infrared continuum flux approaches its maximum, and 
the molecular excitation of high-dipole species should be strongly affected by absorption of 
continuum radiation. Hence any reliable estimation of molecular column densities require 
accurate models for the dust emission. Unfortunately ISO’s lack of angular resolution forces 
us to rely on plausibility arguments in our assumptions about the regions where the different 
lines are formed; some of these are based on general requirements of excitation, and others 
on conclusions from observations of galactic sources. The main goal of this work is thus to 
shed light on the physical and chemical processes in Arp 220, based on detailed model fits 
of its continuum and far-infrared molecular/atomic line absorption and emission spectrum. 
We adopt a distance to Arp 220 of 72 Mpc (projected linear scale of 350 pc/arcsec, Graham 
et al. 1990). In section 2 we present the ISO/LWS observations; in section 3 we discuss 
the line identifications; section 4 is a discussion of the models for the continuum emission; 
section 5 presents the models for the molecular and atomic species; in section 6 we discuss 
the implications of the radiative transfer models, and section 7 summarizes our main results. 


2. Observations 

The full 43-197 /zm spectrum of Arp 220, obtained with the LWS spectrometer (Clegg 
et al. 1996) on board ISO (Kessler et al. 1996) (TDT 2 =27800202), was presented by Fischer 
et al. (1997, 1999). The grating spectral resolution is «0.3 /zm in the 43-93 /zm interval 
(detectors SW1-SW5), and of «0.6 in the 80-197 /zm interval (detectors LW1-LW5), corre- 
sponding to Av > 10 3 kms' 1 . The lines are thus not resolved in velocity. The beam size of 
« 80" ensures that all the continuum and line emission/absorption from Arp 220 (CO size 
< 4", Scoville et al. 1997, hereafter SYB97) lie within the ISO aperture. 

The data we present here was reduced using version 9.1 of the Off Line Processing 
(OLP) Pipeline system which we found to produce a higher signal-to-noise spectrum than 
OLP 10.1. However, we adopted the continuum correction given by OLP version 10.1, which 
typically gives absolute responsivity corrections with uncertainty factors ~ 3 times lower 
than are produced by version 9 (Tim Grundy, private communication). In order to obtain 
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a smooth spectrum throughout the whole LWS range, the density fluxes given by each 
detector were corrected by multiplicative scale factors. Corrections were lower 
than 10% except for detector SWl (43-64 /im), for which the correction was of 
15%. Thus we attribute a conservative uncertainty of 15% to the overall continuum level. 

The LWS spectrum of Arp 220 is presented in Figure 1 together with identifications of 
the most prominent lines. Owing to transient effects, the fluxes of weak lines as observed in 
the forward and reverse scans were found to differ significantly in some wavelength ranges. 
In these cases, if a line appeard close to the upper or lower end of a detector, 
the reverse or forward scan was selected for that line, respectively, to minimize 
the transient effects (Tim Grundy, private communication). Nevertheless, the 
average of both scans was used throughout most of the infrared spectrum. In 
wavelength regions where two detectors’ responses overlap, line fluxes were generally found 
to be consistent. The only exception was the H 2 0 3 22 — 2n line at 90 /zm, which showed 
in LW1 a flux 60% weaker than in SW5. We adopted here the SW5 spectrum, but the flux 
of the above H 2 0 line should be considered highly uncertain. The subtraction of a baseline 
(see Fig. 1) added additional uncertainty to the line fluxes, particularly in cases of broad 
features presumably composed of several lines. With the exception of the H 2 0 3 22 — 2n line, 
we estimate a line flux uncertainty generally lower than 35%. 


3. General results 

OH and H 2 O: The FIR spectrum of Arp 220 is dominated by unresolved OH dou- 
blets (that will be simply referred as lines) and H 2 0 lines in absorption, with 
the exception of the OH 3/2 — 1/2 emission line at 163.3 /zm. Figure 2 shows the level 
diagram of OH, ortho-H 2 0 and para-H 2 0, and indicates the lines detected in Arp 220. Lines 
with very different excitation requirements are observed throughout the spectrum. The OH 
lines n 3/2 J = 9/2 - 7/2 at 65 /zm and Ui/ 2 J = 1/2 — 5/2 at 71 /zm have lower levels at 290 K 
and 415 K above the ground state, respectively, whereas the H 2 0 lines 4 32 — 3 2 i (59 /zm) and 
4 22 — 3u (58 /nn) have lower levels at 305 K and 205 K. Strong absorption is also observed in 
the OH ground state lines at 53, 79, and 119 /zm, as well as in the H 2 0 lowest-lying line at 
179 /zm. This wide range of excitation suggests that several regions with different physical 
conditions are contributing to the observed features (Fischer et al. 1999), and one of the 
goals of this work is to provide a reasonable estimate of the nature of these regions and 
their relative contributions to the spectrum. On the other hand, several of the lines have 
complex shapes, with evidence of shoulders suggestive of weaker secondary lines. In partic- 
ular, the ground state OH 119 /zm line shows a redshifted “shoulder”, which is detected in 
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both the forward and reverse scans, although with somewhat different strengths. It could be 
attributed to the 18 OH II 3 / 2 5/2 — 3/2 line at 120.1 /tm, although contamination by other 
species such as CH + cannot be ruled out. Also, the redshifted “wing” of the H 2 0 2 12 — l 0 i 
line at 179 p m, attributed in Fig. 1 to the H 2 0 2 21 — 2 12 line, could also be contaminated 
by H£ 8 0, CH (see Fig. 7), and H 3 0+. 

CH and NH: The spectrum contains lines from other molecular species: like CH at 149 
pm, NH 3 at 125, 166 and 170 pm and, very interestingly, strong absorptions at 102 and 
153.2 pm that have been identified as NH in Fig 1. Evidence for the latter identifications is 
strengthened because of the presence of weak line-like features at 155.74 pm, and marginally 
at 151.53 pm, which would correspond to the NH 2i — lj and 2i — lo lines, respectively (see 
also Fig. 3). Conceivably, the line absorptions at 102 and 153.2 pm could be severely con- 
taminated by other species, like C 3 , H^O, NH 3 , and even OH + . C 3 has a strong transition 
at 153.3 pm, but its contribution is expected to be minimal due to the lack of detection 
of other adjacent C 3 lines (in particular at 154.86 pm). The absorption at 102 pm may 
be contaminated by the H^O 2 20 — In line (just at 102.0 pm), but since the H^O should 
not be as strong as the corresponding adjacent line of the main isotope H^O, we regard 
this identification as also unlikely. Some contribution of NH 3 lines to the 102 pm feature 
might be expected, but they are somewhat shifted in wavelength (they lie between 101.5 and 
101.7 pm). Both features could be contaminated to some extent by OH + , with strong lines 
at 101.70 and 101.92 pm ( Nj — 3 3 — 2 2 and 34 — 2 3 ) and at 153.0 pm (2 3 — 1 2 ), but the 
strong absorption at 153.2 pm absorption could never be explained by OH + alone. There- 
fore, despite the possible contamination from other molecules, NH is probably responsible 
for most, if not all, of the observed 153.2 pm absorption (see also Fig. 7). Thus, although 
the definitive assignment to NH should await confirmation with higher spectral resolution 
observations of the lowest-lying NH transitions at ~ 10 3 GHz, we conclude that ISO ob- 
servations strongly support its detection in Arp 220, and advance a model of the observed 
absorption (section 5) that can be useful to direct future observational and theoretical stud- 
ies. If confirmed, this detection is the first extragalactic detection of NH, which has been 
previously detected only in galactic diffuse clouds through electronic transitions (Meyer & 
Roth 1991; Crawford & Williams 1997) and, interestingly, toward Sgr B2 in the Galactic 
Center via the same transitions detected in Arp 220 (Cernicharo, Goicoechea, & Caux 2000; 
Goicoechea, Rodriguez- Fernandez, & Cernicharo 2004, hereafter GRC04). 

NH$: The spectrum of Arp 220 around 125 and 170 pm, shown in Figure 4, strongly 
supports the identification of NH 3 . The shape of the 165.7 pm feature indicates transient 
effects, but the line is detected in both the forward and reverse scans. Some H 2 0 lines may 
contribute to the observed 125 and 127 pm features, but they are shifted in wavelength 
relative to the strongest absorption. This is the first extragalactic detection of infrared NH 3 
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lines. With the detection of NH and NH 3 we might expect to detect NH 2 , but there is no 
evidence for its strongest expected lines at 159.5, 117.8 and 104.9 ^ra. 

CO?: An apparent emission line, detected in both the forward and reverse scans, is 
present at 173.7 /xm in Fig. 1. It coincides rather well with the expected position of the CO 
J = 15 — 14 line at 173.63 /xm. This identification cannot be confirmed by the detection 
of other expected CO lines because the CO J — 14-13 line at 186.0 /xm lies at the noisy 
edge of the LW5 detector and the CO J = 16 - 15 line at 162.8 /xm is blended with the OH 
163.3 /xm line. The higher J lines are expected to be too weak to be detectable, given the 
observed strength of the 173.7 /xm feature. 

[O I] and [C II]: The ISO spectrum also shows the [O I] 63.2 /xm line in absorption and 
the [C II] 157.7 /xm line in emission. The [0 I] 145.5 /xm line is not detected. These lines, as 
observed in Arp 220 and other ULIRGs, have been discussed elsewhere (Fischer et al. 1997, 
1999; Luhman et al. 1998, 2003). In section 5 we present a simple model of the Arp 220 
spectrum that may shed some light on the peculiar behavior of these lines in Arp 220. 

In summary, the far-infrared spectrum of Arp 220 shows molecular lines of OH, H 2 0, 
CH, NH, and NH 3 . The atomic lines of [O I] at 63 /xm and [C II] at 158 /xm are also 
detected. Lines of other species, like CO, H 3 0 + , CH + , and OH + , could also contaminate the 
observed features, but our limited spectral resolution prevents the possibility of unambiguous 
detection. Only the [C II] 158 /xm and the OH II 1/2 3/2-1/2 163 /xm lines are clearly observed 
in emission. Lines from ions that would trace H II regions and/or an AGN are absent. 


4. Models for the continuum 

Figure 1 shows that the continuum peaks around ~ 40 - 50 /xm. The bulk of the 
continuum from Arp 220 is emitted by heated dust grains. At 1.3 millimeter wavelengths, 
Sakamoto et al. (1999) showed that the continuum arises almost exclusively from the nuclei, 
with an equivalent size of ~ 0'.'4. The non-thermal contribution at 1.3 mm is expected 
to be < 15% (cf. Fig. 6 of Anantharamaiah et al. 2000). On the other hand, Soifer et 
al. (1999, hereafter S99) have shown that the two nuclei also account for essentially all the 
continuum at 25 /xm. Combining both observations, S99 proposed two alternative scenarios 
to explain the continuum emission of Arp 220 from far-infrared to millimeter wavelengths. 
Our models of the continuum emission are entirely based on these scenarios, which we have 
examined and refined quantitatively on the basis of our ISO 45-200 /xm spectrum. 

Model S\: In the first scenario (hereafter Si), it is assumed that the emission from the 
nuclei is not significantly attenuated at 25 /xm by foreground material. We have simulated 
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the emission from the nuclei as axising from a single nucleus with effective size of 0'.'41. With 
an effective dust temperature of 85 K, and optically thick emission in the submillimeter, the 
requirement that the fluxes at 25 and 1300 /wn arise from the nuclei is fulfilled. However, 
the emission at 60-100 (im , as well as the total luminosity from the galaxy, are then under- 
estimated and a more extended region (hereafter ER) must be invoked to account for the 
remaining flux. We identify this surrounding environment with the extended emission ob- 
served in CO, HCN, and HCO + . In this first scenario, then, most of the Arp 220 luminosity 
is produced in the ER, which has been modeled as a thin disk by SYB97, and as a warped 
disk by Eckart & Downes (2001). Significantly, if this model is correct, a spatially extended 
starburst, responsible for the bulk of the far-infrared luminosity is inferred. 

Our best fit to the continuum using model Si assumptions is presented in Fig. 5a, with 
derived physical parameters listed in Table 1. In all models, uniform densities throughout 
the different components are assumed for simplicity. In Table 1, A t is the wavelength for 
which the nucleus becomes optically thin (r = 1); owing to the high opacities involved, the 
inferred emission is rather insensitive to the spectral index /?. Thus we have given values of 
the physical parameters for /? = 1.5 and P = 2. The dust mass has been derived by assuming 
a mass-opacity coefficient of 12 cm 2 /gr at 200 ^m (Hildebrand 1983). The parameters that 
have been allowed to vary in our models of the ER are the dust temperature T*, the diameter 
d (within the range 500-800 pc), and the dust opacity; /? is fixed to 2 to ensure negligible 
emission at millimeter wavelengths (see S99). 

Model S 2 : In the second scenario (hereafter S 2 ), the emission from the nuclei at 25 /zm is 
assumed to be attenuated by foreground dust with r a b S (24.5 /urn) = 1.2, a value which was 
chosen to be compatible with the silicate absorption observed in S99. The nuclei account for 
the required flux at 24.5 and 1300 /zm with Tj « 106 K, significantly warmer than the ~ 85 
K temperature in Si and, as before, the emission in the submillimeter is optically thick (see 
Table 1). In this scenario as well, however, the flux at 60-100 /zm is again underestimated, 
and an emitting ER must also be involved to account for it. Nevertheless, the luminosity 
from the warm nuclei in S 2 is enough to account for the observed total luminosity from 
Arp 220, so that the ER merely re-radiates the emission from the nuclei and no extended 
starburst is then needed to provide the bulk of the luminosity. 

Figure 5b shows our best fit for S 2 . The unique parameters that have been allowed to 
vary in our models of the ER are the diameter d and the dust opacity; /3 is again fixed to 
2, and the dust temperature T, j (shown in the insert panel of Figure 5b) has been computed 
from the requirement that the heating balance the cooling throughout the source. The 
calculation of Tj is carried on by assuming spherical symmetry, with the nucleus, 
the primary heating source, located at the center of the ER. The ER is divided 
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into a set of spherical shells to account for the variation of T d with the distance to the 
nucleus. Once the T d profile is calculated, the flux contributions of the attenuated 
nucleus and the ER are computed separately, and added up to give the total flux. 
Despite the good fit to the continuum in Fig. 5b, S2 implicitly supposes a lack of spherical 
symmetry (e.g. the nuclear disk by SYB97) or some dumpiness, because the derived 
radial opacity of the ER at 24.5 /xm is 11.3, whereas the adopted opacity of the absorbing 
shell in front of the nuclei is r a bs(24.5pm) = 1.2. Furthermore, we do not rule out 
the possibility that the ER is only partially responsible for the foreground dust 
absorption of the nuclear emission. If the ER were concentrated in a thin disk as 
proposed by SYB97, little dust in the ER would be expected to lie in front of the 
nuclei and significant dust absorption would be attributed to another component, 
“the halo” (see section 5.1). Therefore our results for S 2 , which assume spherical 
shapes for the nucleus and the ER, should be considered only approximate, but 
suggestive. The intrinsic geometry that underlies S 2 departs from spherical symmetry 
and implies that the total luminosity, which coincides with the luminosity of the nuclei, is 
lower than the value inferred in Si, where spherical symmetry and uniformity is assumed for 
each component (Table 1). 

Table 1 shows that /? = 2 yields a dust mass of ~ 10 8 M 0 for the nucleus, and that Si 
gives a mass 1.4 times higher than S2. Since SYB97 infer a dynamical mass of 6 — 8 x 10 9 
M 0 enclosed in the inner 250 pc radius, and this region contains the nuclei and most of the 
inner disk (the ER), S2 with /? = 1.5 is favoured in our models provided that the gas-to-dust 
mass ratio is not lower than the standard value of «100. The consistency between the 
dynamical mass and the mass derived from the dust emission indicates that the 
ISM dominates the dynamics in these inner regions of Arp 220. 

It is worth noting that these models may be applied to the source as a whole, as implicitly 
assumed above, or alternatively to each one of an ensemble of N c smaller clouds of radius 
R c that do not spatially overlap along the line of sight. The value of N c x R 7 C determines the 
absolute scale, and the radial opacity and temperature distribution of each cloud as a function 
of the normalized radial coordinate, R c , determine the continuum shape. Identical results 
are found as long as the above parameters remain constant. Furthermore, both alternatives 
give identical total masses, but differ in the inferred mean density, which scales as y/W c . 
N c = 1 gives the lowest mean density < n(H2 >, which is listed in Table 1. Typical values 
of a fewxlO 4 cm -3 are derived for the nuclei, accounting for the emission from molecules 
with high dipole moment such as HCN, HCO + , HC3N and CN (see also SYB97). For the 
ER we obtain n(H 2 ) < 10 3 cm -3 ; since HCN and HCO + appear to show extended low-level 
emission (Radford et al. 1991b), it is suggested that the actual density is higher than this 
lower limit or that the gas is clumpy. The low density derived for the ER may also 
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be a consequence of the spherical shape attributed to the ER: if the mass we 
derive for the ER were concentrated in a thin disk with full thickness of 32 pc 
(SYB97), the mean H 2 density would be 7 x 10 3 cm -3 . We will adopt the density 
given in Table 1 , n(H 2 ) = 5.3 x 10 2 cm -3 , in the models for molecules and atoms, 
but will also explore the results obtained with a density one order of magnitude 
higher than the quoted value. 

Although both scenarios Si and S 2 reproduce the continuum emission from Arp 220 over 
the 25-1300 pm interval and support the constraints on the nuclear sizes derived from the 
available high angular resolution continuum measurements, the dynamical masses inferred 
from CO millimeter line observations favour S 2 with /? = 1.5 for the nuclei. Moreover, as we 
discuss in sections 5.3 and 6.1, the observed line absorption/emission also favours model S 2 . 
We thus adopt scenario S 2 for the detailed modeling and analysis of the line emission and 
absorption. 


4.1. Extinction 

In both scenarios, the high brightness and compactness of the nuclei imply extreme 
continuum optical depth, corresponding to Ay ~ 10 4 mag. This conclusion is in strong 
contrast with the much more moderate values derived from infrared and radio hydrogen 
recombination lines (Genzel et al. 1998; Anantharamaiah et al. 2000). The high extinction 
derived here is the direct result of the measured 1.3 mm continuum flux from the nuclei, 
210 mJy (Carico et al. 1992; Sakamoto et al. 1999), and the observed upper limit of the 
corresponding source size, ~ 0"4 (Sakamoto et al. 1999). These values imply r d 1 ' 3mm xT d (K) ~ 
40, which shows that even assuming unexpectedly high average dust temperatures (e.g. 
T d = 200 K) and /? = 1 the dust emission is still optically thick even at 200 pm. On the 
other hand, the radio recombination lines observed by Anantharamaiah et al. (2000) are not 
affected by dust obscuration, although their predicted fluxes and the derived extinction may 
be somewhat model dependent. These very different extinction values may be understood 
if we assume that the observed H recombination lines, tracing primarily star formation, 
are formed in the outermost regions of the nuclei, while a buried central energy source, 
responsible for the heating of dust in the innermost regions of the nuclei and a significant 
fraction of the galactic luminosity, is weak in recombination lines. If weak in recombination 
lines, the buried energy source is presumably weak in PAH features and PDR lines as well, 
consistent with the strong [C II] deficit in Arp 220 (see also section 5). Dust-bounded ionized 
regions, in which most of the Lyman continuum from nuclear starbursts or AGN is absorbed 
by dust rather than by gas, may explain these properties of Arp 220, as was proposed by 
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Luhman et al. (2003). We further argue that the Lyman continuum luminosities derived 
from recombination lines do not empirically rule out the possibility that an AGN accounts 
for more than ~ 50% of the bolometric luminosity of Arp 220, because of the high dispersion 
of Iboi/^Lyc values shown by both starburst galaxies and AGN, the range of L L yc values 
derived from different tracers in Arp 220, and the uncertainties in the assumed extinction 
law and the derived extinction (Genzel et al. 1998). Our derived N(H 2 )~ 10 25 cm“ 2 is high 
enough to obscure a source of high 5-10 keV luminosity from one or both nuclei in Arp 220, 
so that a hidden AGN is allowed despite the relatively weak X-ray luminosity observed in Arp 
220 (Clements et al. 2002). Haas et al. (2001) have also argued that a hidden AGN powers 
much of the luminosity of Arp 220 on the basis of the observed submillimeter continuum 
excess relative to the 7.7 \i m PAH flux. 


5. Models for molecules and atoms 
5.1. Comparison with Sgr B2 

The comparison of the spectrum of Arp 220 with that of some well-studied galactic 
sources provides important clues about the regions where the observed lines are formed, 
while emphasizing the unique features that characterize the extragalactic source. In this 
sense, Sgr B2 (component M) is an ideal comparison source, just because it shares common 
observational properties with Arp 220 despite the obvious differences in spatial scale (and 
indeed possibly in nature). Figure 6 shows the continuum-normalized spectra of Sgr B2 
(M) (kindly provided by J.R. Goicoechea) and Arp 220. Sgr B2 harbors newly born OB 
stars, ultracompact H II regions, and hot cores, enshrined in a dusty envelope which is 
heated by shocks and by radiation, and which radiates a high IR luminosity (see GRC04 and 
references therein). The envelope has the highest extinction in the direction of the N- and 
M-condensations, with optically thick emission in the far-infrared up to ~200 /zm, and with 
foreground absorption lines of OH, H2O, CH, and [O I] observed with the ISO/LWS grating 
(GRC04, Fig. 6). Fabry-Perot observations of Sgr B2 have also allowed the detection of other 
molecular species like NH, NH 3 , NH2, HD, H30 + , and C3, as well as high excitation lines of 
OH and H 2 0 (GRC04 and references therein). As shown above (section 4), the extinction 
toward the nuclei of Arp 220 is also very high, although the dust is significantly warmer 
than in the Sgr B2 envelope. The comparison between both sources is meaningful (at least 
as a first approximation, and from the point of view of the radiative transfer) as long as the 
nuclear region of Arp 220 can be considered an ensemble of continuum-thick molecular-rich 
clouds such as Sgr B2. For an ensemble of Sgr B2-like clouds, since both the continuum and 
the line absorption scale with the number of clouds, the continuum-normalized spectrum is 
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the same as that of one individual cloud. This result applies even if the lines in the ensemble 
are broadened relative to the one-cloud emission due to cloud-to-cloud velocity dispersion 
and rotation, provided that the lines remain unresolved with the grating resolution. Thus 
the differences between the two spectra in Fig. 6 reveal real differences in excitation and/or 
column densities. 

The high-excitation lines of OH and H 2 0 are much stronger in Arp 220 than in Sgr 
B2 (M) (see Fig. 6). In particular, the OH II 3 / 2 J = 9/2 — 7/2 65 /zm line, with strong 
absorption in Arp 220, is not detected in the grating spectrum of Sgr B2, and the OH 
n 3/2 J = 7/2 — 5/2 84 /z m line is also much weaker in Sgr B2 (M). This strongly indicates 
the presence in Arp 220 of a high excitation region with relatively high OH column densities. 
Fabry-Perot spectral resolution (Ai> ~35 kms' 1 ) allowed Goicoechea & Cernicharo (2002) 
to detect high excitation OH lines in Sgr B2, and showed that they are pumped through 
absorption of far-infrared photons. In Arp 220 also, these lines appear to be pumped by 
the strong infrared radiation flux in the neighbourhood of the nuclei (section 5.2). Toward 
Sgr B2 (M), Goicoechea & Cernicharo (2002) derived N( OH)« 2 x 10 16 cm“ 2 , and we may 
expect significantly higher column densities toward the nuclei of Arp 220. The peculiarity of 
Arp 220 is also revealed by the relatively strong absorptions in the NH 3 and NH lines. These 
species have also been detected toward Sgr B2 (M) with Fabry-Perot spectral resolution 
(Ceccarelli et al. 2002, GRC04), but Fig. 6 indicates much higher column densities in Arp 
220, at least toward the nuclei. 

GRC04 found that the OH fli/ 2 J = 3/2 — 1/2 163 /zm line, pumped through absorption 
of photons in the OH Hi / 2 — n 3 / 2 J = 3/2 — 3/2 53.3 /zm line, shows emission over a large 
region associated with Sgr B2. In Arp 220 the line is strong, suggesting significant widespread 
emission, i.e. from the ER. It is also worth noting that, although the [C II] line is not detected 
in the grating spectrum of Sgr B2 (M), Fabry-Perot observations allowed its detection (Vastel 
et al. 2002, GRC04), with a flux of « 1.5 x 10” 17 W cm -2 for the component observed in 
emission. This value is, within a factor of 2, similar to strengths in the surrounding region 
where the continuum is however much weaker (GRC04). Therefore, in addition to effects 
of self-absorption and absorption of the continuum by C + in foreground excitation clouds 
(Vastel et al. 2002, GRC04), the low [C II]/FIR ratio at Sgr B2 (M) is due to a strong increase 
in the FIR emission that is not accompanied by a corresponding rise of PDR line emission. 
Dust-bounded ionized regions (Luhman et al. 2003), together with extinction effects in the 
far-infrared, might account for the lack of PDR line emission associated with this additional 
infrared component. This is similar to our hypothesis for Arp 220’s nuclei (section 5.6). The 
[C II] emission in Arp 220 is expected to arise from PDRs in the ER, where the bulk of 
the PAH emission is found (Soifer et al. 2002). As shown in Fig. 5, the ER dominates the 
observed far-infrared emission, but since the main heating sources are, according to model 



-13- 


S 2 , the nuclei, the ER is mainly heated via absorption of infrared rather than UV photons, 
its intrinsic luminosity is relatively low, and hence the [C II] line remains weak (section 5.6). 

Finally, we stress that Arp 220 and Sgr B2 (M) show similar absorptions in the ground- 
state lines of OH and H 2 0, as well as in the CH line at 149 pm. Fabry-Perot observations 
of the low-lying OH and H 2 0 lines toward Sgr B2 (M) (Cemicharo et al. 1997, GRC04), 
indicate that most of these absorptions are produced in diffuse low-excitation clouds located 
along the line of sight to Sgr B2 but not physically associated with it. The similar absorp- 
tions found in Arp 220 strongly suggest that a diffuse medium is also present there. In 
fact we have found that the combination model of the nuclei and the ER that reproduces 
reasonably the high excitation OH and H 2 0 lines, fails to explain the strong absorptions 
in the lowest-lying OH and H 2 0 lines. The presence of an absorbing diffuse component in 
Arp 220 is supported by (z) the detection of CH at 149 pm, because most of 
the corresponding absorption toward Sgr B2, with strikingly similar strength, is 
associated with translucent clouds (Stacey, Lugten, &; Genzel 1987, GRC04) and 
is therefore expected to trace primarily gas in diffuse clouds; {it) the observed 
absorption in the [O I] 63 pm line, which could suggest foreground absorption 
since most of the comparable [O I] absorption observed toward Sgr B2 is also 
produced by foreground gas in diffuse clouds (Baluteau et al. 1997); (iw) the 
low ratio of the CO (2-1) to the 1.3mm continuum emission toward the nuclei 
of Arp 220 and, most important, the low brightness of both the CO (1-0) and 
CO (2-1) lines (SYB97), which suggest line-of-sight blocking of the nuclear CO 
emission by low-excitation gas (Sakamoto et al. 1999). The above points strongly 
suggest the presence of an absorbing component in front of the nuclear region in 
which both the particle density and the infrared radiation density are relatively 
low. This component could also account for significant absorption of the nuclear 
continuum emission in scenario S 2 if the ER were a thin disk (section 4). We 
will refer to this diffuse component of Arp 220 as “the halo” . 


5.2. Outline of the models 

The dust models described in section 4 set up the basis for the molecular calculations. 
These are carried out with the method described in Gonzalez-Alfonso &; Cernicharo (1997, 
1999), which computes the statistical equilibrium populations of a given molecule by assum- 
ing spherical symmetry and line broadening caused by microturbulence and/or radial velocity 
gradients. In the present calculations we have assumed, for simplicity, pure microturbulent 
line broadening, but some tests showed that the inclusion of a radial velocity gradient hardly 
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modified the results. Rotational motion, which is present around the nuclear region of Arp 
220 (e.g. SYB97, Sakamoto et al. 1999; Downes & Solomon 1998), is not included. Never- 
theless, the assumption of spherical symmetry and our neglect of steep velocity gradients 
that may result from cloud-to-cloud velocity dispersion may be considered more critical. 

Our non-local code accounts for radiative trapping in the molecular lines, collisional 
excitation, and excitation through absorption of photons emitted by dust. The dust param- 
eters derived above for Si and S 2 are used in the calculations for molecules. As in the case of 
the continuum models, we have modeled the nucleus and the ER separately. This is required, 
in S 2 , by the relatively low dust opacities of the ER toward the nucleus as compared with the 
radial opacity of the ER (see section 4), and involves inevitably an additional uncertainty. 

Since the dust in the nucleus is very optically thick throughout the ISO wavelength 
range, model results are only sensitive to the molecular column densities in the external- 
most parts of the nucleus. It is therefore assumed that only in the externalmost regions of 
the nucleus, where the infrared lines are formed, the molecular abundances are different from 
zero (see Goicoechea & Cernicharo 2002, for the case of Sgr B2). Since dust and molecules 
are assumed to be coexistent, extinction effects within the nucleus are implicitly taken into 
account; they place important constraints on the molecular abundances. We have adopted 
a molecular shell thickness of 2 x 10 18 cm, which for a mean n(H 2 )= 4.6 x 10 4 cm -3 (in S 2 ) 
corresponds to A v ~ 50 mag and r (50 pm) ss 0.3. For lines around 50 pm, the contribution 
to the absorption by molecules located deeper into the nucleus was checked in some tests to 
be relatively weak, due to both dust and molecular line optical depth effects. Foreground 
extinction in S 2 was also taken into account. In the models for the ERs, which have much 
lower continuum opacities, we have assumed that dust and molecules have uniform abun- 
dance ratio throughout the whole region. The presence of the central nucleus is included in 
the calculation of the statistical equilibrium populations, but ignored in the calculation of 
the emergent fluxes to avoid accounting for it twice. 

Toward the nucleus, absorption of continuum radiation determines the excitation of 
OH and H 2 0; the radiative rates are much higher than the collisional ones even for very 
high pressure molecular gas, such as is found in molecular shocks (n(H 2 )= 5 x 10® 
cm -3 and T k = 300 K; we will refer to these values as “shock conditions”). The 
rate coefficients of Offer, van Hemert, k van Dishoeck (1994) and Green, Maluendes, & 
McLean (1993) were used to check the collisional excitation of OH and H 2 0, respectively. 
If widespread shock conditions were present, only the absorption of the lowest-lying 
lines would be significantly affected. Since we use a halo to match these lines, we cannot 
distinguish between shock and non-shock conditions (i.e., the line ratios are not sensitive to 
n(H 2 ) and T k within plausible values). For simplicity, physical parameters for the halo are 
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derived by assuming non-shock conditions for the nucleus. Concerning the ER, widespread 
shock conditions are not applicable because they would involve strong emission in the H 2 O 
2i2 - loi and 3o3 — 2i2 lines and in most OH lines (hard to cancel by any halo). Our 
molecular data are therefore only sensitive to the radial molecular column density N and 
the microturbulent velocity dispersion a v , so that only these two computational parameters 
are required to define a model of a given component. 

Similar to the dust models, models for molecules may be applied to the source as a whole, 
or alternatively to each one of an ensemble of smaller clouds that do not spatially overlap 
along the line of sight. Besides the scaling relationships pointed out above, the molecular 
column density must remain the same to obtain identical results when varying N c . On the 
other hand, once the number of continuum sources N c is fixed (e.g., N c = 1), the same line 
fluxes are obtained if both N and a v are divided by the same factor / c , and the resulting 
fluxes are then multiplied by /„. The latter reflects the approximate equivalence between 
one absorbing cloud with column density N and velocity dispersion and f c clouds, with 
parameters N/ f c and <?„/ f c , which overlap on the sky but not in the line-of-sight velocity 
space 3 . 

Variations in N and a v have different effects on the line absorptions. If cr v decreases the 
absorptions are weaker for optically thick lines, but due to the increase of line opacities, the 
high-energy levels become more populated and the above weakness is more pronounced for 
the low-lying lines. On the other hand, the increase of N has little effect on very optically 
thick lines (most of them low-lying lines), but a larger effect on those with moderate opacities. 

We have generated a grid of models for the nucleus and the ER of both Si and S 2 , 
by varying the above free parameters N and o v . In each model, the molecular shell is 
divided into a set of sub-shells in order to account for the spatial variation of the excitation 
temperatures of the lines. First we searched for the nucleus+ER combination model that 
best matches the OH and H 2 0 non-ground-state lines, with the same value of o v for both 
species. As pointed out above, the relatively deep absorptions of the OH and H 2 0 ground- 
state lines could not be fitted satisfactorily by any model, and a halo component was added 
to match these lines. The halo was assumed to be a purely absorbing shell; although the 
equilibrium populations were computed in spherical symmetry assuming a size of three times 
that of the ER, limb emission (i.e. emission for impact parameters that do not cross the 
continuum source) was ignored in the calculation of the fluxes emergent from the halo. 

Once the model for OH and H 2 0 is determined, models for the other detected species 


3 In fact the models could be defined in terms of the two independent variables N/cr v and N x / c ; never- 
theless we will use the variables N and a v with f c = 1. 
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are performed, also keeping fixed the value of u v derived above for each component. Figure 
7 compares the derived model spectrum and the observed one, and Table 2 lists the inferred 
parameters. The next sections are devoted to explaining the details of these calculations and 
results. 


5.3. OH and H 2 0 

Our best model fits for the high excitation lines of OH and H 2 0 come from model S 2 , 
and the results presented here will focus on this scenario. The main difference between Si 
and S 2 consists in the higher column densities and/or broader line widths required by Si 
to reproduce the lines, owing to the fact that T& is significantly lower in Si. Models for 
the nucleus with broad line widths (a v > 60 kms' 1 , see below), however, predict strong 
absorptions in some OH and H 2 0 lines, such as the OH ni/ 2 J = 5/2 — 3/2 98 /zm and the 
H 2 0 3 is — 2 02 138 nm lines, which are not observed. Therefore S 2 , which still requires high 
column densities, is a better fit to the data. 

Our best model fits for the high excitation lines involve column densities of 2 — 6 x 10 17 
cm -2 for both OH and H 2 0, and a v = 50 — 30 kms -1 , towards the nucleus. The model for 
the nucleus reproduces nearly the whole absorption in the OH ni/ 2 7/2 — 5/2 and ni/ 2 — 
n 3/2 5/2 — 5/2 lines, most of the OH Fl 3 / 2 9/2 — 7/2, and significant absorption in the other 
lines but by far too weak in the ground-state lines. It also reproduces the full absorptions in 
the H 2 0 4 22 — 3 i 3 , 4 32 - 3 2 i, 330 - 2 21 , 3 3 i - 2 20 , and 4 23 - 3 i2 lines, and significant absorption 
in the 3 2 i — 2i 2 and others. The somewhat low value derived for o v is required to keep the 
absorptions weak in some low excitation lines which are marginally or not detected. Since 
low a v implies low velocity coverage for absorption of the continuum, the column densities 
that are needed to explain the absolute values of the absorptions in high excitation lines 
are relatively high. We stress that this value of a v must be interpreted as a strong lower 
limit on the linewidths that would be observed with high enough spectral resolution, because 
rotation and cloud-to-cloud velocity dispersion would broaden the observed lines. The CO 
kinematic models of Arp 220 by SYB97 found a generic value of a v = 90 kms -1 ; as the 
authors discuss this value should be considered the joint effect of the local linewidth and 
the cloud-to-cloud velocity dispersion over a scale of ~100 pc. Our a v is the local linewidth 
involved in the calculation of opacities and directly related to the molecular excitation, and 
thus the kinematic value of SYB97 must be considered here an upper limit. 

It is worth noting that the above column densities are derived by forcing the OH and 
H 2 0 high excitation lines to arise in the same region. A slight improvement to the fit of 
the H 2 0 lines is obtained with even lower a v , 25 kms -1 , and JV(H 2 0)~ 10 18 cm -2 . This 



-17- 


may indicate that H 2 O and OH toward the nucleus do not arise in the same regions. We 
will adopt in the following the nucleus model with a v — 50 kms" 1 , corresponding to the 
spectrum shown in Fig. 7 and the parameters given in Table 2, and estimate an uncertainty 
of a factor of ~3 on the derived column densities. 

Once the model for the nucleus is matched, the combination nucleus+ER that better 
accounts for the remaining flux in non ground-state lines is searched for. For densities in 
the ER 71 (H 2 ) < 10 4 cm -3 (section 4) and T k = the OH and H 2 0 collisional 
excitation is found to be negligible in comparison with the radiative excitation. 
We have also explored the plausible situation that the OH lines are formed 
within the C + region of PDRs (see section 5.6): assuming T k = 300 K (the 
maximum allowed by the collisional rates of Offer et al. 1994) and n(H 2 ) < 10 4 
cm“ 3 , radiative rates still dominate over collisional rates, and results are found 
indistinguishable from those obtained with lower T k values. Only densities above 
10 5 cm" 3 with Tjt = 300 K would give results significantly different for the ground- 
state transitions of OH. 

The ER mainly accounts for the OH 3/2 - 1/2 163 /im line, which is uniquely 
observed in emission (but predicted in absorption in the nuclei), for more than half of the 
absorption observed in the OH n 3 / 2 7/2 — 5/2 line, for reemission in the 5/2 — 3/2, 
and for significant absorption in the three ground-state OH lines. We estimate for the ER 
iV(OH)~ 2 x 10 17 cm -2 , with a v = 50 kms -1 throughout most of the ER and a v = 90 
kms -1 just around the nucleus. This higher value of a v was required to obtain significant 
reemission in the non-detected OH ni/ 2 5/2 — 3/2 line; since here geometrical effects 
may be important, this result should be considered with caution. For H 2 0 we obtain a 
significantly lower 7V(H 2 0)~ 3 x 10 16 cm" 2 , giving significant absorption in the 3 2i — 3i 2 , 
220 — In, and 2 2 i — lxo lines, and some reemission in the 3o3 — 2 i2 line. Since the halo also 
yields some absorption in these H 2 0 lines but much deeper absorption in the ground-state 
2i2 — loi one (see below), the relative H 2 0 column density in these components is not well 
determined. 

In the halo, the values of N (OH) and o v were determined by fitting the missing ab- 
sorption in the three OH ground-state lines. A value of a v relatively low, 15-20 kms" 1 , and 
JV(OH)~ 2 x 10 16 cm -2 , were found to reasonably fit the fluxes of the 79 and 119 /im lines, 
though the flux of the 53 fim line is somewhat underestimated (Fig. 7). Higher values of 
cr v would predict too much absorption in the already saturated 119 /im line. For H 2 0 we 
find 7V(H 2 0)~ 1-5 x 10 16 cm" 2 . Owing to the strong radiation field from the nucleus and 
ER, H 2 0 in the halo is still significantly excited, thus yielding also some absorption in the 
220 — In, 2 2 i — lio and 3 q3 ~ 2i 2 lines. 
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Despite the generally satisfactory fit obtained for the OH and H 2 0 lines, the shoulder 
of the OH 119 /z m line, presumably produced by the 18 OH n 3 / 2 5/2 — 3/2 line, is not 
reproduced. The models assume 16 OH/ l8 OH=500; however, the isotopic abundance ratio in 
Arp 220 required to reasonably fit the 120 /xm shoulder is as low as 16 OH/ 18 OH~50. This 
value cannot be ruled out if 14 N is efficiently converted into 18 0 in nuclear processing of 
high mass stars, and then efficiently ejected to the interstellar medium through stellar winds 
and/or supernovae (Henkel & Mauersberger 1993). Also, there is compelling evidence for 
isotopic ratios of 150-200 in starburst regions of nearby galaxies (Henkel & Mauersberger 
1993). Due to the low spectral resolution of the spectrum, and the possibility that the feature 
is contaminated by other species, we do not attempt to place useful constraints on this ratio. 
Nevertheless, we conclude that a very low 16 OH/ 18 OH abundance ratio likely applies to Arp 
220, perhaps indicating an advanced stage starburst (Henkel & Mauersberger 1993). 

The model of Fig. 7 predicts an absorbing flux of 1.8 x 10 -12 erg s -1 cm -2 for 
the II 1/2 — n 3 /2 5/2— 3/2 line at 34.6 /zm, in reasonable agreement with the observed 
flux of 2.1 x 10 -12 erg s' 1 cm -2 reported by Skinner et al. (1997). Concerning 
the OH-megamaser emission in the 1.667 GHz line, our models predict inversion 
in the line, but with low gain and by far too weak to account for the observed 
emission. Proper models for the maser emission require a much finer spatial grid 
than that used in the present calculations, and we have not attempted to account 
for it. Even so, some remarks can be given on this score. It is suspected that the 
OH-megamaser emission is radiatively pumped through absorption of photons in 
the 34.6 /zm and 53.3 /zm lines, followed by radiative cascade to lower levels. Since 
the lower level of the 34.6 /zm transition is the ground-state n 3 / 2 J = 3/2 level, we 
have found that about 65% of the modeled absorption is predicted to occur in the 
foreground halo, rather than in the nuclei. As a consequence, the pump efficiency 
that would be required in our model to explain the 1.667 GHz OH-megamaser 
line via the 34.6 /zm line alone is relatively high, ~ 2%. Nevertheless, given the 
uncertainty we have in the derived nuclear OH column density, the possibility 
that some of the maser emission arises in more inner regions, and the expected 
additional contribution to the pumping by the 53.3 /zm line, we conclude that 
our models are roughly consistent with the OH-megamaser excitation scheme 
discussed by Skinner et al. (1997). 


5.4. CH 


Since CH is close to the Hund’s coupling case (b) limit in its 2 n ground 
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state (e.g. Brown &: Evenson 1983), we denote its levels through (iV, J), where 
N is the case (b) rotational quantum number and J = N ± We assume that the 
(2,3/2)-(l,l/2) CH line we observe at 149.2 /z m arises in the halo, based on the results 
obtained toward the Galactic Center (see GRC04). We derive 7V(CH)« 2 x 10 15 cm“ 2 by 
fitting the feature. Unlike the case of Sgr B2, however, the submillimeter emission from 
the nuclear region of Arp 220 is strong enough to populate significantly the (1,3/2) 
level via absorption of photons in the (1,3/2)— (1,1/2) line at 560 /zm, so that our 
models for the halo predict some contribution by CH (2,5/2)-(l,3/2) to the absorption 
feature at 181 /zm (Fig. 7). The latter is uncertain, however, because the feature at 181 /zm 
could be contaminated by H 3 0 + (GRC04) and/or by stronger absorption of H^O, whose 
abundance relative to H^O could be enhanced relative to the assumed value of 1 /500. On 
the other hand, there is a wing-like feature at 118.5 /im, observed in both the forward and 
reverse scans, which could be caused by the doublet CH (3,7/2)-(2,5/2). If so, and since the 
excitation of this line requires a relatively strong radiation field, there would be CH in the 
nuclei that would account for about 1/3 of the absorption at 149.2 /z m, and the CH column 
density in the halo would be 2/3 of the quoted value. 


5.5. NH and NH 3 

In contrast with OH and CH, most of the column density of NH and NH 3 we 
model is contained within the nuclei. Assuming that the 153.2 ^m absorption feature is 
caused entirely by NH, we have obtained N{ NH)~ 10 16 cm' 2 toward the nucleus. The model 
reproduces the marginal absorption feature at 76.8 /zm, attributable to the Nj = 4j — 3 j* 
lines. For the above column density, the 22 — li and 2 3 — 1 2 lines at 153.2 /xm become 
saturated and the associated feature is not completely reproduced, so that we have added 
an additional halo component with N( NH)« 2 x 10 15 cm' 2 (Table 2). Nevertheless, this 
model still underestimates the absorption at 102 /zm, strongly suggesting the contribution 
from other species (section 3). 

Within a given AT— ladder, the excitation of the NH 3 non-metastable levels is determined 
by absorption of far-infrared continuum photons, while the metastable levels in K = 2, 3, . . . 
are populated through collisions (see e.g. Ceccarelli et al. 2002, for the case of Sgr B2). 
In the model of Fig. 7 we have assumed an average density of n(H 2 )= 4.6 x 10 4 cm' 3 and 
T k = 100 K, but we have checked that the model results are insensitive to the adopted 
I* because of the blending of lines from different K — ladders to each spectral feature in 
our spectrum (Fig. 1). We obtain N(NH 3 )~ 3 x 10 16 cm -2 toward the nucleus to fit the 
absorptions at 125 and 127 /xm (Fig. 7). Besides the absorption in the lines showed in 
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Fig. 4, the model predicts significant absorption at wlOO /xm and «101.6 /xm, caused by 
(J, K) = (5, K) - (4, K) lines, which is still insufficient to account for the observed 
absorption around 101.6 /im. The model also fails to explain the strong absorption 
at 166 /jm, and therefore an halo component with A(NH 3 )~ 4 x 10 15 cm -2 has been added 
to the global model of Fig. 7. Nevertheless, the halo components of NH and NH 3 should be 
considered uncertain, because variations in the background continuum associated with each 
component could in principle account for the missing flux in the lines. 

Finally, we have explored the possibility that the NH 2 radical contributes to the spec- 
trum at some wavelengths. The expected strongest absorption from NH 2 is found at ~159.5 
/zm, caused by the strongest components of the 3i 3 — 2 02 ortho line (the hyperfine structure 
was neglected in these calculations, but the split of the levels due to the unpaired electronic 
spin of 1/2 was taken into account). At this wavelength, a marginal absorption feature 
may be attributed to NH 2 , and is approximately fit with a model for the nucleus where 
AT(NH 2 )~ 10 15 cm -2 . We have included this model of NH 2 in Fig. 7 to show that the 
expected absorption in other lines, like the 3 22 — 2 U one at 105 /mi, do not conflict with the 
observations, and we conclude that AT(NH 2 )< 2 x 10 15 cm -2 . 


5.6. C II and O I 

A crucial test of our model is whether it can reproduce the [C II] 157.7 //m emission and 
[O I] 63.3 /mi absorption lines. Among ULIRGs, Arp 220 shows one of the most extreme [C 
II] deficits (F cn /F FIR « 2xl0' 4 , Luhman et al. 2003). The [C II] line is formed within A v <2 
mag from the surfaces of PDRs (e.g. Wolfire, Tielens, & Hollenbach 1990), where the UV 
field from nearby high mass stars, or from the average galactic field has not been significantly 
attenuated. In this region, photodissociation maintains most of the gas in atomic or singly 
ionized form, but some radicals, like OH and NH, find their maximum abundances there 
(Sternberg & Dalgarno 1995). In particular, OH is expected to be an excellent molecular 
tracer of PDRs’ surfaces, given that its abundance is rather low in UV-shielded quiescent 
molecular clouds. Its abundance relative to H nuclei within the C + region of dense PDRs 
is expected to approach the value of ~ 3 x 10 -6 (Sternberg & Dalgarno 1995). In fact, 
Goicoechea & Cernicharo (2002) have found an abundance of « 2 x 10 -6 in Sgr B2, and 
it could be as high as 5 x 10~ 6 around the galactic center (Genzel et al. 1985). We have 
estimated A(OH)~ 1 - 3 x 10 -6 toward the nucleus of Arp 220. 

On the above grounds, and adopting a gas phase carbon abundance of 1.4 x 10 -4 (Savage 
& Sembach 1996), we assumed N(C + )/A(OH)=100 in Arp 220 (Table 2), and we computed 
the expected [C II] line emission by assuming excitation through collisions with atomic 
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H (Tielens & Hollenbach 1985). The collisional rates were taken from Launay & Roueff 
(1977a). The H densities were assumed to be twice the H 2 densities derived from 
the continuum models, i.e. n(H) = 9.2 xlO 4 cm -3 in the nucleus and n(H) = 1.06x10 s 
cm -3 in the ER. Since the critical density is 3 x 10 3 cm' 3 (Kaufman et al. 1999), results for 
the nucleus are not critically dependent on the assumed density. They are also insensitive 
to the assumed temperature as long as it is higher than ~100 K (Wolfire et al. 1990). Thus 
our results are only sensitive to the assumed density in the ER, and to the assumed C + 
abundance. 

The result of this calculation has also been added to the overall model of Fig. 7 (7* = 500 
K has been assumed). The model overestimates the observed [C II] emission in only 
24%. The contribution from the nucleus is only 13% of that from the ER because of the 
low volume of the nuclear emitting region. The expected line flux from the nucleus could be 
lower if absorption of the underlying continuum by low-excitation C + in the halo, ignored 
in this model, occurs as observed in Sgr B2 (M) (GRC04). The bulk of the line emission 
arises from the ER. If the density of the ER were one order of magnitude higher 
than assumed (section 4), the modeled [C II] emission would be a factor of « 2 
stronger than in Fig. 7. The situation toward Arp 220 resembles what is found in Sgr B2 
(GRC04), where the line is emitted mainly from an extended region around condensations 
N and M, while the strong FIR source itself is not associated with corresponding observable 
[C II] line emission. 

The [O I] 63.3 /xm line has been modeled by assuming iV(0 I) /AT (OH) =250 (Table 2); 
although the oxygen abundance is expected to be twice that of C + in the atomic region, it 
is expected that further atomic oxygen exists deeper into the clouds (Sternberg &; Dalgarno 
1995). For this reason, a wide range of excitation temperatures is expected for the [O I] line. 
We have just fitted a single “effective” kinetic temperature and assumed also collisions with 
atomic H. The same densities as assumed above for C + excitation are used for O I, and the 
collisional rates are taken from Launay & Roueff (1977b). 

Our calculations show that absorption in the [O I] 63.3 fim line is obtained, both toward 
the nucleus and the ER, with an effective T* = 160 K, but the line is still too weak to 
account for the observed feature. Given the high O I abundances that are expected in 
diffuse clouds (Baluteau et al. 1997), we have added to the model a halo component with 
iV(0 I)= 3 x 10 18 cm -2 (Table 2). This model accounts for the observed absorption at 63.3 
Hm (Fig. 7). 

On the basis of the low extinction derived from infrared and radio H recombination lines 
and the high optical depth derived from our dust models (section 4.1), and the assumption 
that [C II] emission is expected to suffer extinction similar to that of the recombination 
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line emission, we favor the explanation proposed by Luhman et al. (2003, see also section 
5.1) that non-PDR far-infrared emission is responsible for the [C II] deficit in Arp 220. 
Nevertheless, our continuum models by themselves (section 4) cannot rule out the possible 
role of far-infrared extinction on the measured line fluxes: the derived high nuclear far- 
infrared opacities indicate that only the externalmost regions of the nuclei, where the OH 
and H 2 0 lines are formed, and the ER, are able to contribute to the [C II] line emission. Our 
model fits cannot discern whether the intrinsic nuclear [C II] emission is negligible or rather 
is obscured by dust. In either case, we have shown that the [C II] line is well reproduced by 
assuming that C + and OH are coexistent. 


6. Discussion 

6.1. The Extended Region (“ER”) 

Our models support the widespread presence of PDRs in the ER. Both the high OH and 
C II column densities indicate that the UV field from newly formed stars have a profound 
effect on the chemistry in the ER. Significant contributions from shocks can be neglected, as 
was pointed out in section 5.2. The H 2 0-to-OH abundance ratio is significantly lower than 
1, probably indicating enhanced H 2 0 photodissociation. As pointed out in section 4, 
the density is uncertain, in the range n(H 2 ) = 5 x 10 2 — 7 x 10 3 cm -3 . 

Our models also indicate that the bulk of the [C II] line emission arises in the ER. It 
is therefore likely that star formation is responsible for this emission. This result is strongly 
supported by the observations of Soifer et al. (2002), who found that the PAH emission is 
also spatially extended. Assuming a “normal” [C II]/FIR ratio, i.e. F( 0 n] /E fir = 5 x 10 -3 
(Stacey et al. 1991), the expected intrinsic FIR emission from the ER is ~ 3 x 10 10 L 0 , i.e. 
~3% of the total galactic luminosity. This estimate supports the scenario S 2 that has been 
used to model the line emission. According to our models, the total FIR luminosity (due to 
absorption and re-emission of nuclear infrared radiation) from the ER is much higher than 
the intrinsic (PDR) emission. Our estimate of the intrinsic ER luminosity is somewhat lower 
than that derived by Soifer et al. (2002) and Spoon et al. (2004), who found a lower limit of 
~ 7 x 10 10 L 0 , but confirms qualitatively their results. The starburst luminosity in the ER 
seems to be similar to that of moderately bright infrared galaxies, like NGC 253 (Radovich, 
Kahanpaa, & Lemke 2001). 

The warm mass in the ER derived from our OH models is ~ 10 9 M 0 , assuming 
X(OH)=10 -6 relative to H nuclei and the warm H 2 arising from the same volume as the 
OH. This estimate of mass may be considered an upper limit. From the H 2 rotational lines 
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detected by ISO, Sturm et al. (1996) estimated ~ 3 x 10 9 M© of warm gas. Given the 
uncertain correction by extinction, it is likely that an important fraction of the H 2 emission 
arises in the ER. 


6.2. The halo 

The models also indicate the presence of a halo. Despite the uncertainties in the column 
densities of this component, they are typical of those found toward our galactic nucleus: 
jV(H 2 0, OH)~ 2 x 10 16 cm -2 have been also derived in the diffuse medium toward Sgr B2 
(Cernicharo et al. 1997; Neufeld et al. 2003; Goicoechea &; Cernicharo 2002), and the derived 
7V(CH)« 2 x 10 15 cm -2 is also similar to that found toward Sgr B2 by GRC04 and Stacey et 
al. (1987). The column densities in the halo derived for NH and NH 3 are uncertain because 
they are based on single lines, but they could also exist in a population of molecular clouds 
located far away from the nuclear region. Assuming X(H 2 0)= 10 -6 cm -2 (Neufeld 
et al. 2003) in the halo, iV(H 2 ) ~ 2 x 10 22 cm -2 is obtained and the associated 
continuum opacity is T 200Mm ~ 8 x 10 -3 . If we further assume an spectral index 
/3 = 2, the derived dust opacity at 25 /x m is ~ 0.5, comparable to the value of 
1.2 in scenario S 2 (section 4). Therefore, significant absorption of the nuclear 
continuum emission is attributable to the halo. 


6.3. The nuclei 

The similarity of the OH column densities in the ER and the nucleus may suggest 
that, at least to some extent, we are observing the same widespread OH component, and 
that OH is more excited toward the nucleus because of the underlying stronger infrared 
continuum in that direction. On this ground a PDR origin of the observed OH would be 
favoured. However, the inferred high H 2 0 column density would be difficult to explain in this 
context. Although H 2 0 is expected to form efficiently in UV-shielded regions of dense PDRs 
(Sternberg & Dalgarno 1995), with total column densities similar to those of OH, there seems 
to be no clear correlation between A' (OH) and Af (H 2 0). In fact, the H 2 0 column density in 
the ER is significantly lower than that of OH. Also, detection of NH 3 and, above all, of NH, 
seem to indicate an additional nuclear component. 

OH could also arise in C-shocks (Watson et al. 1985), but it is unlikely that they dom- 
inate the OH absorptions because Af(H 2 0) would then be at least one order of magnitude 
higher than jV(OH). However, an interesting possibility is that those C-shocks, or alterna- 
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tively hot core regions, are combined with PDRs, i.e., H 2 O produced there is subject to a 
strong UV field and thus to photodissociation. This process may be enhanced if there is hot 
gas emitting X-rays, like in supernova remnants where the OH abundance is expected to be 
> 10 -6 (Wardle 1999), or from a nuclear AGN. In the nuclei of Arp 220, several compact 18 
cm continuum sources indicate the presence of high luminous supernovae (Smith et al. 1998); 
however, the extended and external OH required to explain the infrared data suggests a more 
widespread component. The diffuse OH megamaser emission found in Arp 220 should be re- 
lated to it (Lonsdale et al. 1998). Soft-extended and hard-compact X-ray emission, detected 
around and from the nuclei (Clements et al. 2002), could be responsible for photodissociation 
of H 2 O produced in shocks and hot cores, thus enhancing the OH abundance. In particular, 
hot cores are expected to exist widely in the nuclei, given the high dust temperatures and 
densities found there. The presence of J-shocks, where 7V(OH) is expected to be higher than 
7V(H 2 0) except for high enough preshock densities (Neufeld & Dalgarno 1989), cannot be 
disregarded. 

The high column densities obtained for NH and NH 3 seem to indicate that standard 
gas-phase PDR chemistry alone is not able to explain the full molecular data in Arp 220. 
Sternberg & Dalgarno (1995) predicted for a PDR an NH column density more than two 
orders of magnitude lower than that of OH, whereas we estimate /V(OH)//V(NH)~20 in Arp 
220. The enhancement of NH relative to OH may be more than one order of magnitude. In 
Sgr B2, GRC04 also found a somewhat high iV(NH) relative to OH, i.e. iV(OH)/iV(NH)~30- 
100. It is interesting that the high NH abundance in diffuse clouds is a factor of ~ 40 higher 
than predicted by gas-phase chemical models (Meyer & Roth 1991). The latter has been 
used by Wagenblast et al. (1993) and Crawford & Williams (1997) to argue for grain-surface 
production of NH. In principle, this process could also help to enhance the NH abundance 
in Arp 220, because the dust in the nuclei has been found to be warm so that grain mantles 
could efficiently evaporate, and also because in an enviroment with enhanced cosmic rays, 
as expected from a starburst and the consequent high rates of supernovae, or an AGN, the 
release of grain mantles to the gas phase by means of sputtering should be also enhanced. 
However, hydrogenation of NH should in principle continue until saturation, because indeed 
we observe NH 3 with a column density twice of NH. The issue that now arises is that, if 
hydrogenation generally completes, then a very low NH-to-NH 3 abundance ratio would be 
expected, and if it does not complete like in the models of Wagenblast et al. (1993), the 
scenario fails because of the low relative abundance found for NH 2 . In Sgr B2, for example, 
NH:NH 2 :NH 3 =1:10.T00 (GRC04), and NH 2 is found to be fairly abundant (van Dishoeck et 
al. 1993). 

One possible solution for the low NH 2 abundance is that the already invoked enhance- 
ment of cosmic rays deeply affects the ion-molecule gas-phase chemistry. Federman, Weber, 
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& Lambert (1996) have shown that the high NH abundance found in some diffuse galactic 
enviroments could be explained through cosmic ray ionization of atomic nitrogen, followed 
by hydrogen abstraction reactions that form NH+ and dissociative recombination that yields 
NH. The slight endothermicity of N + +H 2 — >NH + +H (Millar, Farquhar, &; Willacy 1997) is 
not a problem here, given the high dust temperatures in the nuclei. If in Arp 220 the cosmic- 
ray ionization of N is enhanced, the above scheme could give rise to high NH abundances. 
Furthermore, both NH 2 and NH 3 would be much less abundant if photodissociation is im- 
portant in those regions. NH 3 would be formed primarily in grain mantles through nearly 
complete hydrogenation, thus again keeping the NH 2 abundance low, and released to the 
gas phase in widespread hot core regions relatively shielded from UV fields. H 2 0 could also 
follow this last process. 

A chemistry deeply influenced by ion-neutral reactions have been also invoked by Aalto 
et al. (2002) to explain the high emission from HNC relative to HCN in Arp 220 and other 
luminous infrared galaxies. Furthermore, Aalto et al. (2002) also found strong subthermal 
CN emission, indicative of gas at moderate densities and irradiated by UV fields (Rodriguez- 
Franco, Martin-Pintado, & Fuente 1998). In the simplest scenario, the molecular content 
of Arp 220 may thus be interpreted in terms of hot cores submitted to strong UV and X-ray 
fields, where enhanced evaporation of grain mantles and ion-molecule chemistry induced by 
cosmic ray ionization are also deeply affecting the relative molecular abundances. 


7. Summary 

We have analyzed the ISO/LWS spectrum of Arp 220 using radiative transfer models 
applied to both the continuum and line emission. Our main results are: 

1. The continuum emission from 25 to 1300 pm is well reproduced with a two-component 
model: (a) the nuclei, with effective size of 0'.'4 and dust temperature of 106 K, which accounts 
for essentially the whole flux at 25 pm and at millimeter and submillimeter wavelengths, and 
(b) an extended region (ER), whose effective size is 2" and which dominates the continuum 
emission from 60 to 250 pm. 

2. The extinction toward the nuclei is very high ( Ay ~ 10 4 mag); the dust in the ER is 
heated through absorption of radiation emanating from the nuclei. 

3. The spectrum of Arp 220 shows molecular lines of OH, H 2 0, CH, NH, and NH 3 , as well 
as the atomic [O I] 63 pm line in absorption and the [C II] 158 pm line in emission. The 
outermost regions of the nuclei, along with the ER, are traced by the lines observed in the 
far-infrared. 

4. The high excitation lines of OH and H 2 0 are pumped through absorption of photons 
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emitted by dust. Column densities of iV(OH)= 2 - 6 x 10 17 cm -2 and iV(H 2 0)= 2 — 10 x 10 17 
cm -2 are derived toward the nuclei. In the ER, iV(OH)~ 2x 10 17 cm -2 and iV(H 2 0)~ 3x 10 16 
cm -2 . We found it necessary to invoke a third component, or halo, to match the low-lying 
lines of OH and H 2 0; this halo has column densities that are similar to those found toward 
the Galactic Center (A(OH, H 2 0)~ 1.5 x 10 16 cm' 2 ). 

5. The CH line detected in the far-infrared spectrum of Arp 220 is assumed to arise from the 
halo, and the inferred column density is iV(CH)~ 2 x 10 15 cm -2 . This value is also similar 
to that found toward the Galactic Center. 

6. Models for NH and NH3 indicate high column densities toward the nuclei, 7V(NH)~ 
1.5 x 10 16 cm* 2 and 7V(NHa)~ 3 x 10 16 cm* 2 . The upper limit found for the column density 
of NH 2 is much lower, iV(NH 2 )< 2 x 10 15 cm* 2 . 

7. The [C II] 158 pm line strength is approximately reproduced by assuming that C + is 
100 times more abundant than OH. Our models predict that the line arises mainly from 
the ER, and that non-PDR far-infrared emission, with possible extinction effects, is mostly 
responsible for the observed [C II] deficit in Arp 220. The [O I] 63 pm line is also matched 
with an abundance of 250 relative to OH and absorption toward the nuclei, the ER, and the 
halo. 

8. PDR molecular chemistry plays a key role in the ER and contributes to the elevated 
OH abundance at the expense of H 2 0. Toward the nuclei, however, important contributions 
from hot cores, and possibly from shocks, is most plausible. The nitrogen chemistry, and in 
particular the high NH abundance, seems to be strongly influenced by ion-neutral reactions 
triggered by cosmic ray ionization. 


We are grateful to J.R. Goicoechea for providing the ISO/LWS spectrum of Sgr B2 and 
for fruitful discussions. E. G-A thanks Spanish SEEU for funding support under project 
PR2003-0057, and the Harvard-Smithsonian Center for Astrophysics for its hospitality. This 
work has been partially supported by NASA Grant NAG5-10659, the NASA LTSA program 
and the Office of Naval Research. 
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Fig. 2. — Energy level diagrams of OH and H 2 O (ortho and para). Rotational levels with 
energies above the ground state up to 620 K are shown; the lines detected in Arp 220 are 
indicated with arrows and their wavelengths are in gm. OH A-doubling is ignored because 
the A-doublets are not resolved with the ISO grating resolution. 




Fig. 3. — Line assignments in the vecinity of the 153.2 /im feature. The grey line shows the 
adopted continuum level. 



Fig. 4. — NH 3 lines around 125, 127, 166, and 170 n m. The grey line shows the adopted 
continuum level. The 125 and 127 /J,m features could be partially contaminated by the 
labelled H 2 0 lines. 
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Fig. 5. — Fits of the continuum emission from Arp 220 in the 25-1300 nm range for scenarios 
(a) Si and (b) S 2 . Solid line shows ISO-LWS spectrum of Arp 220, and filled triangles show 
the 60 and 100 fim IRAS fluxes for comparison. Filled squares show the fluxes measured by 
Soifer et al. (1999, 24.5 fim); Eales, Wynn-Williams, & Duncan (1989, 450 /xm); Rigopoulou, 
Lawrence, & Rowan-Robinson (1996, 350, 800, and 1100 /xm); and Sakamoto et al. (1999, 
1300 jxm). Dotted lines indicate the computed contributions from the nuclei and the ER, 
whereas the solid grey line show the expected total flux. The insert panel shows the calculated 
dust temperatures in function of the radial angular distance for S 2 (see text). 
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Continuum-normalized spectra of Sgr B2 (M) and Arp 220. The main carri< 
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Table 1. Physical parameters derived from the continuum emission 
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Table 2. Column densities derived toward Arp 220* 


Species 

Nucleus 

ER 

Halo 

OH 

2.0 x 10 17 

2.0 x 10 17 

1.8 x 10 16 

H 2 0 

2.0 x 10 17 

3.0 x 10 16 

1.5 x 10 16 

CH 

- 

- 

2.0 x 10 15 

NH 

1.3 x 10 16 

- 

1.8 x 10 15 

nh 3 

2.8 x 10 16 

- 

4.0 x 10 15 

C+ 

2.0 x 10 19 

2.0 x 10 19 

- 

0 

5.0 x 10 19 

5.0 x 10 19 

2.7 x 10 18 


“Units are cm 
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ABSTRACT 

We report on the analysis of the first complete far-infrared spectrum (43- 
197pm) of the Seyfert 2 galaxy NGC 1068 as observed with the Long Wavelength 
Spectrometer (LWS) onboard the Infrared Space Observatory (ISO). In addition 
to the 7 expected ionic fine structure emission lines, the OH rotational lines 
at 79, 119 and 163pm were all detected in emission, which is unique among 
galaxies with full LWS spectra, where the 119pm line, when detected, is always 
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in absorption. The observed line intensities were modelled together with ISO 
Short Wavelength Spectrometer (SWS) and optical and ultraviolet line intensities 
from the literature, considering two independent emission components: the AGN 
component and the starburst component in the circumnuclear ring of ~ 3 kpc 
in size. Using the UV to mid-IR emission line spectrum to constrain the nuclear 
ionizing continuum, we have confirmed previous results: a canonical power-law 
ionizing spectrum is a poorer fit than one with a deep absorption trough, while the 
presence of a big blue bump is ruled out. Based on the instantaneous starburst 
age of 5 Myr constrained by the Br 7 equivalent width in the starburst ring, 
and starburst synthesis models of the mid- and far-infrared fine-structure line 
emission, a low ionization parameter (U=10~ 3,5 ) and low densities (n=100 cm -3 ) 
are derived. Combining the AGN and starburst components, we succeeded in 
modeling the overall UV to far-IR atomic spectrum of NGC 1068, reproducing 
the line fluxes to within a factor 2.0 on average with a standard deviation of 
1.3, and the overall continuum as the sum of the contribution of the thermal 
dust emission in the ionized and neutral components. The OH 119 pm emission 
indicates that the line is collisionally excited, and arises in a warm and dense 
region. The OH emission has been modeled using spherically symmetric, non- 
local, non-LTE radiative transfer models. The models indicate that the bulk of 
the emission arises from the nuclear region, although some extended contribution 
from the starburst is not ruled out. The OH abundance in the nuclear region is 
expected to be ~ 10~ 5 , characteristic of X-ray dominated regions. 

Subject headings: galaxies: individual (NGC 1068) — galaxies: active — galaxies: 
nuclei — galaxies: Seyfert — galaxies: starburst — infrared: galaxies. 


1. INTRODUCTION 

NGC 1068 is known as the archetypical Seyfert type 2 galaxy. It is nearby, luminous 
(Lir = 2 x lO n L0 Bland-Hawthorn et al. 1997), and it has been extensively observed and 
studied in detail from X-rays to radio wavelengths. With a measured redshift of z=0.0038 
(Huchra et al. 1999) (corresponding to a distance of D=15.2 Mpc for H 0 =75 km s -1 Mpc -1 ), 
it provides a scale of only ~ 74 pc/". A central nuclear star cluster has an extent of ~ 


'ISO is an ESA project with instruments funded by ESA Member States (especially the PI countries: 
Fiance, Germany, the Netherlands and the United Kingdom) and with the participation of ISAS and NASA. 
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0.6" (Thatte et al. 1997) and a 2.3 kpc stellar bar observed in the near-IR (Scoville et al. 
1988; Thronson et al. 1989) is surrounded by a circumnuclear starburst ring. Telesco et al. 
(1984) found that the infrared emission in NGC 1068 was due to both the Seyfert nucleus 
(which dominates the 10pm emission) and to the star forming regions in the bright ~ 3kpc 
circumnuclear ring (which emits most of the luminosity at A > 30 pm). A Brj imaging 
study (Davies, Sugai k Ward 1998) showed a similar morphology and indicated that a short 
burst of star formation occurred throughout the circumnuclear ring of 15-16" in radius within 
the last 4-40 Myr. CO interferometer observations revealed molecular gas very close to the 
nucleus (~0.2") suggesting the presence of ~ 10 8 M© within the central 25pc (Schinnerer et 
al. 2000). Recent high resolution H 2 line emission mapping indicates the presence of two 
main nuclear emission knots with a velocity difference of 140 km/s, which, if interpreted as 
quasi-keplerian, would imply a central enclosed mass of 10 8 M© (Alloin et al. 2001). 

In this article, we present the first complete far-infrared spectrum from 43 to 197pm 
showing both atomic and molecular emission lines (§2). We model the composite UV- to 
far-IR atomic emission line and continuum spectrum, from our data and the literature, using 
photoionization models of both the active nucleus and the starburst component (§3). We 
also model the mid- to far-IR continuum emission using a radiative transfer code and gray 
body functions for the neutral molecular components (§4). Moreover, two different non- 
local, non-LTE radiative transfer codes have been used to model the OH lines (§5). Our 
conclusions are then given in §6. 


2. OBSERVATIONS 

NGC 1068 was observed with the Long Wavelength Spectrometer (LWS) (Clegg et al. 
1996) on board the Infrared Space Observatory (ISO) (Kessler et al. 1996), as part of the 
Guaranteed Time Programme of the LWS instrument team. The full low resolution spectrum 
(43-197pm) of NGC 1068 was collected during orbit 605 (July 13, 1997). Two on-source 
full scans (15,730 seconds of total integration time) and two off-source (6’ N) scans of the 
[CII] 158pm line (3,390 seconds of total integration time) were obtained. On- and off-source 
scans had the same integration time per spectral step. Because of the design of the LWS 
spectrometer, simultaneously with the 158pm data, a short spectral scan of equal sensitivity 
to the on-source spectrum was obtained at sparsely spaced wavelengths across the LWS 
range. 

The LWS beam is roughly independent of wavelength and equal to about 80 arcsec. 
The spectra were calibrated using Uranus, resulting in an absolute accuracy better than 
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30% (Swinyard et al. 1996). The data analysis has been done with ISAP 2 , starting from 
the auto-analysis results processed through the LWS Version 7-8 pipeline (July 1998). To 
be confident that newer versions of the pipeline and calibration files did not yield different 
results, we have compared our data with the results obtained using pipeline 10.1 (November 
2001) and we did not find significant differences in the line fluxes or the continuum. 

All the full grating scans taken on the on-source position and the two sets of data on 
the off-source position were separately co-added. No signal was detected in the off-source 
coadds. The emission line fluxes were measured with ISAP, which fits polynomials to the 
local continuum and Gaussian profiles to the lines. In all cases the observed line widths 
were consistent with the instrumental resolution of the grating, which was typically 1500 
km/sec. The integrated line fluxes measured independently from data taken in the two scan 
directions agreed very well, to within 10%. The on source LWS spectrum that resulted from 
stitching the ten LWS channels together using small multiplicative corrections in order to 
match the overlapping regions of each channel with its neighbors is shown in Fig. 1. LWS 
spectra of sources that are very extended within the instrument beam or that peak off center 
are typically affected by channel fringing in the continuum baseline (Swinyard et al. 1998). 
Fortunately, these spurious ripples are hardly noticeable in our LWS spectrum, presumably 
because the far-IR continuum is centrally concentrated towards the center of the LWS 80" 
beam. 

Besides the LWS observations, we also use the SWS observations presented by Lutz 
et al. (2000), to extend the wavelength and ionization-level coverage. Table 1 presents all 
the ISO line flux measurements including those from the SWS with their respective aperture 
sizes. 


3. THE FINE STRUCTURE LINES 

To be able to better constrain the modeling of the line emission of NGC 1068, we have 
combined our far-infrared fine structure line measurements (Table 1) with ultraviolet, optical 
and infrared spectroscopic data from the literature (Kriss et al. 1992; Marconi et al. 1996; 
Thompson 1996; Lutz et al. 2000). The complete emission line spectrum of NGC 1068 from 


2 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are Centre d’Etude Spatiale des Rayonnements (Prance), 
Institute d’Astrophysique Spatiale (France), Infrared Processing and Analysis Center (United States), Max- 
Planck-Insitut fur Extraterrestrische Physisk (Germany), Rutherford Appleton Laboratories United King- 
dom) and the Space Research Organization, Netherlands. 
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the ultraviolet to the far-IR includes several low-ionization lines that are primarily produced 
outside the narrow line region (NLR) of the active nucleus, as well as intermediate ionization 
lines that originate from both starburst and AGN emission. For this reason, we find that 
no single model satisfactorily explains all the observed emission lines. We identify two main 
components: 

- an AGN component (the NLR), exciting the high ionization lines and contributing 
little to the low-to-intermediate ionization lines; 

- a starburst component in the circumnuclear ring of the galaxy (e.g. Davies, Sugai & 
Ward 1998) that produces the low ionization and neutral forbidden lines and some of 
the emission in the intermediate ionization lines. This component should also produce 
emission associated with photo-dissociation regions (PDRs) (e.g. Kaufman et al. 1999), 
at the interface with the interstellar medium of the galaxy. 

In this section, we will examine separately the two components that produce the total 
fine structure emission line spectrum of NGC 1068, namely the AGN and the starburst, for 
which we propose two different computations, and we add together these components to 
reproduce the overall observed spectrum from the UV to the far-IR in §3.3. 


3.1. Modeling the AGN 

The first photoionzation model predictions of the mid to far-infrared emission line spec- 
tra of the Narrow Line Regions (NLR) of active galaxies were presented by Spinoglio & 
Malkan (1992), well before the ISO observations could be collected. Alexander et al. (2000) 
used the observed high ionization emission lines to model the obscured ionizing AGN con- 
tinuum of NGC 1068 and found that the best-fit spectral energy distribution (SED) has a 
deep trough at 4 Rydbergs, which is consistent with an intrinsic “big blue bump” that is 
partially obscured by ~ 6 x 10 19 cm -2 of neutral hydrogen interior to the NLR. Following 
their results, we have simulated their models, although using a different photoionization 
code, CLOUDY (Version 96, Ferland et al. 1998; Ferland 2000), and then we have varied 
the shape of the ionizing continuum to include the ionizing continuum derived in Pier et al. 
(1994). Our goal was to test if the Alexander et al. (2000) results were unique and to fit the 
remaining emission by a starburst component, and thereby to derive a composite model of 
the complete emission line spectrum of NGC 1068. 

Specifically, we explore three plausible AGN SEDs. Model A assumes the best fit ionizing 
spectrum derived by Alexander et al. (2000), i.e. with a deep trough at 4 Rydberg (log / = 
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-27.4, -29.0, -27.4, -28.2 at 2, 4, 8 and 16 Ryd, respectively). An intrinsic nuclear spectrum of 
NGC 1068 has also been inferred by Pier et al. (1994). Model B assumes the original ionizing 
spectrum derived from Pier et al. (1994). Model C assumes an SED with a Big Blue Bump 
superposed on the Pier et al. (1994) ionizing continuum (log / = -25.8, -25.8, -25.8, -27.4 at 
2, 4, 8 and 16 Ryd, respectively) as expected for the thermal emission of an accretion disk 
around a central black hole. These three AGN ionizing continua are plotted in Fig. 1. For 
each of models A, B, and C, we have used two component models with the same parameters 
as in Alexander et al. (2000): component 1 has a constant hydrogen density of 10 4 cm“ 3 4 , 
an ionization parameter U=0.1, a covering factor c=0.45, a filling factor of 6.5 x 10 -3 with 
a radial dependence of the form r -2 , and extends from ~ 21 to ~ 109 pc from the center; 
component 2 has a density of 2 x 10 3 cm -3 , an ionization parameter U=0.01, a covering 
factor of c=0.29, a filling factor of 6.5 x 10 -4 without any radial dependence, and extends 
from ~ 153 to ~ 362 pc from the center. We have also assumed the “low oxygen” abundances 
adopted by Alexander et al. (2000), in order to be able to compare our results with theirs 3 . 
Because the grain physics has been updated in the most recent version of CLOUDY (Version 
96), we have included the presence of grains in the models, using “Orion-type” grains 4 . The 
inclusion of grains also allows us to compute the thermal dust continuum emission from the 
ionized components (see §4). 

The inner and outer radii of the emission regions of the two components, 21, 109, 153 
and 362 pc, correspond to angular distances of about 0.26, 1.4, 1.9 and 4.5 ", respectively. 
Table 2 reports the predicted line fluxes of the three AGN models, A, B, and C, together 
with the observed line fluxes: the line fluxes are given for each of the two components 1 and 
2, which are treated as independent, and the total flux for each model is simply the sum of 
the fluxes of the two components. 

We can see from Table 2 that only the AGN A and B models, and not the AGN C 
model, reproduce most of the observed high ionization line fluxes. The low and intermediate 
ionization lines, are expected to have partial or full contributions from starburst and PDR 
components (see §3.2). This first result rules out the presence of a “big blue bump” in 
the ionizing continuum of NGC 1068. To be able to compare the modeled ultraviolet and 
optical lines with the observations, we also listed in Table 2 and 3 their dereddened fluxes, 


3 The adopted gas phase chemical abundances in logarithmic form are: H: 0.00, He: -1.00, Li: -8.69, Be: 
-10.58, B: -9.21, C:-3.43, N:-3.96, O : -3.57, F : -7.52, Ne: -3.96, Na: -5.67, Mg: -4.43, Al: -5.53, Si: -4.46, 
P : -6.49, S : -4.79, Cl: -6.72, Ar: -5.60, K : -6.88, Ca: -5.64, Sc: -8.83, Ti: -6.98, V : -8.00, Cr: -6.33, Mn: 
-6.54, Fe: -4.40, Co: -7.08, Ni: -5.75, Cu: -7.79, Zn: -7.40 

4 The abundances of the grain chemical composition, in logarithmic form, are: C: -3.6259, O: -3.9526, Mg: 
-4.5547, Si: -4.5547, Fe: -4.5547 



This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article 
may difTer from this preprint. 

Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/'428495". 

-7- 

assuming two values for the extinction: E b-v = 0.4 mag (Malkan & Oke 1983) and E b-v 
— 0.2 mag (Marconi et al. 1996). We find that the AGN B model overpredicts several 
of the intermediate ionization lines, such as [SIV]10.5/xm, [NeIII]15.6/xm and [SIII]18.7/xm, 
and this discrepancy increases when adding the starburst component because these lines 
are also copiously produced by that component (see next section). On the other hand, the 
[NeII]12.8/xm emission is underpredicted so much so that even with the inclusion of the 
starburst component it cannot be reproduced with this model. As we discuss further in 
§3.3, a composite AGN/starburst model using AGN model A reproduces the [NeII]12.8/xm 
emission better than the other composite models. 


3.2. Modeling the starburst ring 

NGC 1068 is known to emit strong starburst emission from the ring-like structure at 
a radial distance of 15 — 16" from the nucleus (total size of ~ 3 kpc), traced for example 
by the Br 7 emission (Davies, Sugai & Ward 1998). Mid-IR line imaging observations 
of NGC 1068 have been published by Le Floc’h et al. (2001) based on ISOCAM CVF 
observations. They presented an image of the 7.7 /mi PAH feature that shows constant 
surface brightness above the 4th coutour near the nucleus. This suggests that star formation 
is occuring in the direction of the nucleus so that nuclear spectra will include some emission 
from star formation. In the case of the SWS observations that we are modeling (reported by 
Lutz et al. (2000)), three apertures were used at different wavelengths with the two largest 
also including portions of the brighter starburst ring (see Table 1). To estimate how much 
of the starburst emission is contained in the different apertures used in the observations, 
we have used a continuum subtracted image in the 6.2 /xm feature produced by C. Dudley 
(private communication) using the same ISOCAM CVF data set examined by Le Floc’h et 
al. (2001). The 6.2 /xm feature is more isolated than the 7.7 /xm feature, which is blended 
with the 8.7 /xm feature and the silicate absorption feature, but the image compares well with 
the published 7.7 /xm image though we have zeroed out the residuals in a 3x 12 arcsec 2 region 
centered on the nucleus. Based on this image, the SWS 14x20, 14x27 and 20x33 arcsec 2 
slits contain 13, 23 and 46% of the 6.2 /xm flux contained in the LWS beam respectively, 
without correction for the neglected region of poor residuals (oriented at 45° to our synthetic 
SWS slits). Since PAH features are thought to be a good tracer of PDRs and their associated 
startbursts, we adopt these percentages in our model predictions of SWS line strengths in the 
starburst models presented in this section. In fitting our starburst models to the observations, 
we have computed the line fluxes at earth of each centrally illuminated emitting cloud and 
then determined the number of clouds needed to best fit the observed line fluxes. 
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We have chosen the starburst synthesis modeling program Starburst99 (Leitherer et 
al. 1999), to produce input ionizing spectral energy distributions (SEDs) for the CLOUDY 
photoionization code. Models were followed to temperatures down to 50 K to include the 
PDR components. We compared the predictions of an instantaneous star formation law with 
those of a continuous star formation law. For both types of models we adopted an age of 5 
Myr, a Salpeter IMF (a=2.35), a lower cut-off mass of 1 Mq, an upper cut-off mass of 100 
M©, solar abundances (Z=0.020) and nebular emission included. These ionizing SEDs are 
shown in Fig. 2 with a total mass of M = 10 6 M© for the instantaneous model and a star 
formation rate of 1 M© yr -1 for the continuous model. These particular ionizing continuum 
shapes were selected because they are consistent with the Br 7 equivalent width observed 
by Davies, Sugai & Ward (1998) in the starburst ring. We have estimated that the Br 
7 equivalent width in each of the individual regions of the map of Davies, Sugai Sz Ward 
(1998) is in the range 110-180 A. According to the Leitherer et al. (1999) models (see their 
figures 89 and 90), for a value of log(W(Br 7, A)) > 2 only instantaneous models with ages 
less than ~ 6x 10 6 yrs are allowed. 

We report in Table 3 the line fluxes predicted for six different centrally illuminated 
starburst models, chosing the above instantaneous star formation SED as the input ioniz- 
ing continuum and using CLOUDY with densities of n h = 10, 100, 1000 cm -3 , ionization 
parameters of Log U = -2.5, -3.5 and an inner cloud radius of 50 pc. As a function of the 
adopted density, we then determined the following numbers of emitting clouds needed to fit 
the observations: 33000, 3300, and 330 clouds for the three values of the density, respective- 
ly. The fluxes reported in Table 3 are therefore the total starburst line fluxes at earth and 
together with the nuclear line fluxes of Table 2 can be compared with the observations. 

We have also rim models with the continuous star formation law presented above, but 
we do not list their results in Table 3, because the differences in the line flux predictions 
compared with the instantaneous models are insignificant, compared with the effects of 
density and ionization parameter, as can be seen from Table 3. This result is not surprising 
because the ionizing continua of the two starburst models are quite similar in shape and 
the total number of clouds is a free parameter. We have also tried continuous starburst 
models with much longer ages (10, 20 and 100 x 10 6 years) but, because the shape of the 
ionizing continuum again does not change significantly, the resulting emission line spectrum 
was indistinguishable from that one derived from the models with an age of 5 x 10 6 years. 

In all models the abundances were those typical of HII regions 5 and grains of “Orion- 


5 The adopted gas phase chemical abundances in logarithmic form are: H: 0.00, He: -1.02, Li: -10.27, Be: 
-20.00, B: -10.05, C: -3.52, N: -4.15, O: -3.40, F: -20.00, Ne: -4.22, Na: -6.52, Mg: -5.52, Al: -6.70, Si: -5.40, 
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type” 6 are included. The integration was allowed to run until the temperature of the gas in 
the cloud cooled to T=50 K in order to include the photodissociation regions present at the 
interfaces of HII regions and molecular clouds. 

It is clear from Table 3 that the models with the higher ionization parameter (log U = 
-2.5) can easily be ruled out, because their emission in many intermediate ionization lines is 
far too high (see e.g. [01 V] 26pm, [OIII] 5 1,88pm, [NIII]57pm). Among the models with the 
lower ionization parameter (log U = -3.5), we can exclude model SBR F, with density n# = 
1000 cm -3 , because it underestimates many far-IR lines which are not strongly emitted by 
the active nucleus (namely: [SiII]35pm, [NIII]57pm, [OIII]88pm, [Nil] 122pm, [CII] 158pm) 
while the low density model (SBR B, with n H = 10 cm -3 ) overpredicts the [CII] 158pm line 
by a factor of 2 relative to the other far-IR lines and does not reproduce the [OIII] doublet 
ratio. Finally, the intermediate density model (SBR D, with n h = 100 cm -3 ) gives the best 
fit to the observed lines, taking into account that the AGN component must be added to 
reproduce the total flux as is shown in §3.3. 

We estimate the average PDR parameters using the models of Kaufman et al. (1999) 
and the contour plots in Luhman et al. (2003), the measured [C II] 158 and [O I] 145 pm line 
fluxes (but not the [O I] 63 pm line flux which may be affected by absorption and/or shocks), 
and the FIR flux integrated over the LWS spectrum, which we find to be 1.3 x 10 -8 ergs 
cm -2 sec -1 . Here we assume that the [C II] line emerges predominantly from PDRs due 
to the strong starburst, rather than the diffuse ionized medium. With this assumption, the 
average PDR gas density and UV radiation field are n# 2 ~ 1000 and Go ~ 300 respectively. 

We note that if instead we assume that the [C II] line flux is dominated by the diffuse ionized 
medium, using the correction factor estimated by Malhotra et al. (2001), we obtain a similar 
gas density n h 2 ~ 1500 but a significantly higher interstellar radiation field Go ~ 1500. For 
both cases, the parameters derived are in the range of those of the normal galaxies in the 
Malhotra et al. (2001) sample, consistent with the assumption that most of the FIR flux 
originates in the starburst ring. 


P: -6.80, S: -5.00, Cl: -7.00, Ar: -5.52, K: -7.96, Ca: -7.70, Sc:-20.00, Ti: -9.24, V:-10.00, Cr: -8.00, Mn: 
-7.64, Fe: -5.52, Co:-20.00, Ni: -7.00, Cu: -8.82, Zn: -7.70 

6 The abundances of the grain chemical composition in logarithmic form are: C : -3.3249, O : -3.6516, Mg: 
-4.2537, Si: -4.2537, Fe: -4.2537. 
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3.3. Adding the two components 

Summing the line intensities of each one of the two components, the composite spectrum 
of NGC 1068 can be derived and compared with the observed one. We have chosen three 
combinations to compute the composite models, each one with a different AGN model, while 
we adopted the starbust model with n h — 100 cm -3 and Log U = -3.5: 1) the first one (that 
we name CM1, for Composite Model 1) with the AGN ionizing continuum as suggested by 
Alexander et al. (2000) (model AGN A); 2) the second (CM2) with the original Pier et al. 
(1994) (model AGN B); 3) the third (CM3) with the hypothetical bump (model AGN C). 
The results of these three composite models are given in Table 4, compared to the observed 
and dereddened values, assuming the two choices for the extinction (see §3.1). We also show 
the results of the three composite models in a graphical way in Fig. 3, where the modeled to 
the observed flux ratio is given for each line for the case of an extinction of Eb-v=0.2 mag. 

A simple y 2 test of the three models resulted in a reduced y 2 of 11.6, 17.1 and 177 for the 
three models CM1, CM2 and CM3, respectively, for an extinction of Eb-v=0.4 mag, while 
these values become 23, 46 and 325 for Eb-v= 0.2 mag. Thus, of the models explored, CM1 
with Eb_v=0.4 mag provides our best fit to the observations. We note that model CM1 
reproduces the line fluxes to within a factor of 2.0 on average, with a standard deviation of 
1.3. 


4. THERMAL CONTINUUM SPECTRUM 

To model the total mid- and far-infrared thermal dust continuum of NGC 1068 we have 
used different computations for each component. The modelled emission is shown for the 
individual components and combined in Fig. 4. The thermal dust emission from the AGN 
narrow line regions and the starburst regions in the ring have been computed using the 
CLOUDY photoionization models described in the previous sections. Specifically, the UV 
continuum reprocessed by the dust present in both NLR components 1 and 2 of model AGN 
A has been diluted by the same covering factors that affect the line emission (c=0.45 and 0.29 
for components 1 and 2, respectively). Similarly we computed the continuum emission from 
our best starburst model (SBR D). While we find that the continuum produced in this way for 
the AGN is consistent with the observed mid-infrared energy distribution (Fig. 4), accounting 
for about half of the observed emission at mid-infrared wavelengths, the emission from dust 
associated with the starburst ionized and photodissociated regions, although similar in shape 
to the observed continuum, produces only ~ 20% of the far-IR continuum. We have indeed 
performed a search in parameter space by varying the age of the starburst (from 4 to 6 x 
10 6 yrs), the gas density (from 10 to 1000 cm -3 ) and the radius of the emitting clouds (from 
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25 to 100 pc), but no starburst model that could reproduce simultaneously the observed line 
and far-infrared continuum emission was found. 

If the CLOUDY models correctly reflect conditions in both the ionized and photodis- 
sociated gas, these results may imply that the bulk of the starburst far-infrared continuum 
arises from dust that is mixed with neutral gas not directly associated with the ionized or 
photodissociated gas. However, because photoionization codes such as CLOUDY have not 
been used in the past to model the dust continuum, and because we may not have fully 
searched parameter space, we are hesitant to over-interpret these results until more detailed 
comparison with galactic ionized and photoionized regions are carried out. We have therefore 
described this starburst thermal dust component, following Spinoglio, Andreani k Malkan 
(2002), in terms of a gray body function with a temperature of T = 34 K, and a colder gray 
body component at T = 20 K (see Fig. 4). These are gray body functions with a steep (/ 3 = 

2) dust emissivity law. Assuming a spherical shell with radius of 1.5 kpc and thickness of 0.3 
kpc, the inferred average H 2 molecular density associated with the 34 K component is 3.7 
cm -3 . The total mass is 1.2 x 10 9 M G and 2 x 10 9 M G for the 34 K and 20 K components, 
respectively. These estimates are in reasonable agreement with the ~ 4 x 10 9 M© derived 
for the molecular ring from CO emission (Planesas et al. 1991) . 

As pointed out above, the dust associated with the ionized NLR components 1 & 2 
of model AGN A does not quite account for the total mid-infrared continuum (Fig. 4). 
Therefore, we have assumed that the missing mid-infrared arises from the neutral-nuclear 
component , which has been observed in a variety of molecular lines (e.g. Tacconi et al. 1994; 
Usero et al. 2004). For this neutral component we have assumed a total gas 4-dust mass of 
M — 2 x 10 7 M 0 (Heifer & Blitz 1995), and modelled the expected continuum using a non- 
local, spherically symmetric, radiative transfer code (Gonzalez- Alfonso k Cernicharo 1997, 
1999). The molecular nuclear emission has been resolved into a circumnuclear disk or ring 
(Schinnerer et al. 2000), which here is roughly modelled as a dusty spherical envelope with 
inner and outer radii of 3 and 200 pc, respectively. We assume an AGN luminosity of 3.7 x 10 10 
L©, which is a factor of ~ 3 lower than the total AGN luminosity, simulating that most of 
the AGN radiation escapes through the poles of the molecular disk, and/or is absorbed in 
the NLR, and so is not able to heat the molecular gas (e.g. Cameron et al. 1993). The dust 
envelope is divided into a set of spherical shells where the dust temperature is computed 
assuming that heating and cooling are equal. We used a standard silicate/ amorphous carbon 
mixture with optical constants given by Draine (1985) and Preibisch et al. (1993). The 
density profile was assumed to be oc with (3 regarded as a free parameter. 

The resulting mid-infrared emission, obtained with (3 — 1, is shown in Fig. 4. The 
nuclear molecular component has an averaged H 2 density at the inner radius of < n/(H 2 ) > = 
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500 cm -3 , and a column density of N( H 2 ) = 2 x 10 22 cm -2 . Once the emission from this 
neutral component is summed up with the emission predicted for the NLR components 1 &; 

2, a good fit to the observed emission for A > 9 /xm is obtained. For A < 9 /xm, the mid- 
infrared emission is underestimated, suggesting the presence of a hot component, probably 
very close to the central AGN, which is not included in our models. 


5. THE OH LINES 
5.1. General remarks 

In NGC 1068, we detect three of the OH rotational lines, all in emission. As shown in 
the energy level diagram of Fig. 5, two of them are fundamental lines, connecting the ground 
state 2 n 3/2 3/2 level with the 2 n 3 / 2 5/2 (the in-ladder 119 /xm line) and with the 2 Hi/ 2 1/2 
level (the cross ladder 79 /xm line). The third line is the lowest transition of the 2 IIi/ 2 ladder: 
the 163/xm line between the J=3/2 and J=l/2 levels. The detected line fluxes are given in 
Table 1. The fact that these three lines are all in emission is in striking contrast with the 
OH lines observed in other bright infrared galaxies, such as Arp 220 (Fischer et al. 1999; 
Gonzalez- Alfonso et al. 2004), Mrk 231 (Harvey et al. 1999), NGC 253 (Bradford et al. 1999), 
and M 82 (Colbert et al. 1999), in which the 119 /xm fundamental is in absorption. The 79/xm 
line is sometimes seen in emission and sometimes in absorption; the 163/xm line is always 
seen in emission. In addition to the detections, the ISO-LWS and SWS observations provide 
upper limits on fluxes of the other four lines that arise between the lowest six rotational 
levels. The LWS spectra in the vicinity of the detected lines (and of one of the upper limits), 
are shown in detail in Fig. 6 (histograms). In this section we discuss the physical conditions 
necessary to excite these lines, their probable location within NGC 1068, and detailed model 
fits to their fluxes. A comparison between the observed and modeled line fluxes is given in 
Table 5 and shown in Fig. 6. 


5.2. The excitation mechanism of the OH lines 

The unique OH emission line spectrum of NGC 1068 can provide a powerful way to help 
discriminate between the properties of the molecular clouds in NGC 1068 and the clouds 
in other galaxies in which OH has been observed. Before describing our detailed radiative 
transfer calculations, it is instructive to discuss some conclusions that are model-independent. 
The emission in the OH n 3 / 2 5/2 — 3/2 line at 119 /xm cannot be explained by absorption of 
far-infrared photons followed by cascade down to the upper n 3 / 2 5/2 level of the transition. 
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Rather, we argue that collisional excitation dominates. Figure 5 shows the energy level 
diagram of OH. There are only two possible paths to excite the 119 /tm line via absorption 
of far-infrared photons: via the 35 /tm and/or the 53 /tm ground-state lines. Excitation by 
either of these routes has other observable consequences. In the case of simple radiative 
cascading, the Einstein - A coefficients of the lines involved in the corresponding cascades are 
such that if the 35 pm absorption path were responsible for the observed 119pm line flux, 
then the OH Iiy 2 5/2 — 3/ 2 line at 98.7 pm would be approximately 5 times stronger than 
the 119 pm line, while the 98.7 pm line is not detected. Hence this possibility is ruled out. 
Similarly, if absorption in the 53pm line were responsible for the observed 119pm line flux, 
then the 163pm line would be about 5 times stronger than the 119 pm line, which it is 
not. We can therefore conclude from the constraints provided by the other far-infrared OH 
lines that the 119pm emission line is not the result of radiative absorption and cascading. 
The implication is that OH excitation through collisions is more important in NGC 1068 
than in the other observed galaxies and therefore that the gas responsible for the observed 
emission in the 119pm line resides predominantly in relatively dense and warm enviroments 
in comparison with these other sources. 

The other two observed emission lines, unlike the 119 pm line, need not be collisionally 
dominated. In the case of the Hi/ 2 3/2— 1/2 163pm line, the most likely excitation mechanism 
is absorption of photons emitted by dust in the 53 and 35 pm lines followed by radiative 
cascade. The upper level of this transition is 270 K above the ground state (Fig. 5), so 
that excitation through collisions is expected to be ineffective in this line. The excitation 
mechanism of the - n 3 / 2 1/2 — 3/2 79 pm line could be a mixure of collisional and 
radiative pumping. The upper level of this transition is 182 K above the ground state, so 
that a warm and dense region could, at least partially, excite the line through collisions. 
Nevertheless, the line could be also excited through the same infrared pumping mechanism 
that results in the observed 163 pm line emission. 

In conclusion, the 119 pm line is collisionally excited, whereas absorption of photons 
emitted by dust in the 53 and 35 pm lines probably dominates the excitation of the 163 pm 
line. The 79 pm OH line may in principle be excited through both mechanisms. 


5.3. Constraints on the spatial origin of the 119 pm OH line 

In NGC 1068 two regions with very different physical conditions can account for the 
observed OH emission as discussed above: the compact nuclear region, and the ring and bar 
where the starburst is taking place. A warm and dense region is required to account for 
the observed 119 pm line emission, given that the line is collisionally excited, so the warm 
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and dense neutral nuclear region around the AGN should be considered a good candidate, 
despite its small size (~ 5"; e.g. Planesas et al. 1991; Schinnerer et al. 2000), for the following 
reasons: 


(z) It is warm: there are ~ 10 3 M 0 of hot H 2 (~ 2000 K) distributed over ~ 5" (Blietz et 
al. 1994) that is thought to be UV or X-ray heated (Rotaciuc et al. 1991). Both PDRs 
and XDRs can produce a range of temperatures as high as a few times 10 3 K (Kaufman 
et al. 1999; Sternberg & Dalgarno 1995; Maloney et al. 1996). From the CO (4-3) 
to (1-0) line intensity ratio, Tacconi et al. (1994) derive ~ 80 K for the bulk of the 
molecular gas, with a mass of ~ 3 x 10 7 M 0 enclosed in within the central 4" (Heifer & 
Blitz 1995). Lutz et al. (2000) have reported the detection of pure H 2 rotational lines 
within the ISO-SWS aperture, and estimated ~ 2.5 x 10 7 M 0 at ~ 200 K, but these 
lines may also arise, at least partially, from the inner regions of the 3 kpc starburst 
ring. 

(zz) The molecular clouds within the nuclear region are dense, although there is some 
dispersion in the values derived by several authors based on HCN emission: Tacconi et 
al. (1994) derived an H 2 density of ~ 10 5 cm” 3 , whereas subsequent observations and 
analysis by Heifer & Blitz (1995) yielded a density of ~ 4 x 10 6 cm” 3 . An intermediate 
density of 5 x 10 5 cm 3 from HCN and CS, and lower for other tracers, has been 
recently derived by Usero et al. (2004). 

(Hi) The OH abundance is expected to reach high values in regions exposed to strong 
incident UV fields (PDRs, Sternberg &; Dalgarno 1995), and in particular in X-ray 
dominated regions (XDRs, Lepp & Dalgarno 1996). The remarkable chemistry found 
by Usero et al. (2004) in the circumnuclear disk of NGC 1068 is indicative of an overall 
XDR and suggests a high OH abundance in the nuclear region. 


Given that the OH 119/zm line is collisionally excited, the possibility that the line might 
arise from the nuclear region can be checked by computing the amount of warm gas required 
to account for the observed emission: 


M w ( M 0 ) = 1.6 x 10 7 x 


10- 5 

v 

'5 x 10 5 cm -3 ' 

X 

4.3 x 10 u cm 3 s 1 

X(OH) 


n(H 2 ) J 

< C/u > 


( 1 ) 


where A (OH) is the OH abundance relative to H 2 , and < c\ u > is the collisional rate for 
excitation from the ground n 3 / 2 3/2 level to the n 3 / 2 5/2 one. Equation 1 assumes that, 
although the line could be optically thick, it is effectively optically thin, and makes use of 
the observed flux of 1.2 x 10” 12 erg s” 1 cm” 2 . The reference value for the collisional rate, 
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< Qti >= 4.3 x 10~ n cm 3 s~\ corresponds to gas at 80 K (Offer et al. 1994); it decreases 
by a factor of « 2.7 for gas at 50 K and increases by a factor of 3 for gas at 200 K. 

The reference OH abundance we use in this estimate, 10~ 5 , is the result of two separate 
studies: first, calculations of molecular abundances by Lepp &; Dalgarno (1996) have shown 
that the OH abundance in XDRs is expected to be about two orders of magnitude higher 
than the abundance of HCN and HCO + . The authors in fact suggested the possibility that 
the high HCN/CO ratio observed in the nuclear region of NGC 1068 could be a consequence 
of enhanced X-ray ionization. Second, the possibility of a chemistry dominated by X-rays 
has found support from observations by Usero et al. (2004), who derive abundance ratios 
of HCN, HCO" 1- , and CN in general agreement with predictions for XDRs. Since the HCN 
abundance derived by Usero et al. (2004) is ~ 10” 7 , X(OH) in XDRs could attain values 
as high as 10 -5 . On the other hand, the density of 5 x 10 5 cm" 3 derived by Usero et al. 
(2004) has been adopted as the reference value in eq. 1. Finally, the mass of molecular gas 
derived from the emission of several molecular tracers is expected to be at least ~ 2x 10 7 M 0 
(Heifer & Blitz 1995), which is similar to the value required in eq. 1. From these estimates 
we conclude that, if the OH abundance is as high as ~ 10 -5 (i.e. if the predictions for XDRs 
are applicable to the nuclear region of NGC 1068), the bulk of the OH 119 /mi line could 
arise there. This possibility would naturally explain why NGC 1068 is unique in its 119 fim 
line emission among galaxies with full LWS spectra. 

Finally we ask whether the OH 119 fim line could arise from an even more compact 
region, i.e., from a torus with a spatial scale of 1 pc surrounding the central AGN. According 
to typical parameters given by Krolik &: Lepp (1989), a torus is expected to be hot (~ 10 3 
K), could have densities of 10 7 cm” 3 , and therefore a mass of ~ 10 5 M 0 . Also, the OH 
abundance is expected to be very high, 5 x 10~ 5 — 10 -4 . Eq. 1 shows that the relatively 
low mass of the torus (about 2 orders of magnitude lower than the entire nucleus) could be 
compensated by the the higher density, OH abundance, and temperature expected there, so 
that this possibility cannot be neglected. 

The reference values for the nuclear abundance and density adopted in eq. 1 are rather 
uncertain (and possibly extreme). The continuum models of §4 indicate that the mass 
associated with the 34 K dust component, which is identified with the starburst ring, is 
1.2 x 10 9 M 0 . If ~ 5% of this mass corresponds to warm molecular gas rich in OH, the 
amount of extended warm gas is ~ 6 x 10 7 M 0 . According to eq. 1, the OH emission at 119 
fim can then also be explained as arising in the ring if the associated PDRs, with assumed 
OH abundance of 2 x 10” 6 (Sternberg & Dalgarno 1995; Goicoechea &; Cernicharo 2002; 
Gonzalez-Alfonso et al. 2004), have densities of a few x 10 5 cm” 3 . Since Papadopoulos &; 
Seaquist (1999b) found that most of the extended molecular gas resides in dense, compact 
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clouds, this scenario seems also possible. However, the continuum from the starburst at 119 
pm is strong, so that one expects that eq. 1 is in this case underestimating M w , and the 
quoted physical parameters, A (OH) and n(H 2 ), are lower limits. The effect of dust emission 
is discussed in detail below. 

In conclusion, a definitive answer to the issue of the spatial origin of the OH 119 pm 
emission cannot be inferred from only the flux observed in the 119 pm line. Nevertheless, 
useful constraints on this subject are given: the line could be either explained as arising from 
the nucleus, with a required OH abundance ~ 10 — 5 , or from the extended ring, with OH 
abundance > 2 x 10 -6 and density > a few x 10 5 cm -3 . Nevertheless, the radiative transfer 
models described below, which take into account the effect of the continuum emission and 
the excitation of the 79 and 163 pm lines, point towards a nuclear origin of the OH emission. 


5.4. Outline of the models 

Analysis of the OH 79 and 163 pm lines requires the use of detailed radiative transfer 
calculations since, as pointed out above, the emission in these lines is expected to be strongly 
influenced by absorption of far-infrared continuum photons. We therefore proceeded to model 
the OH lines with two different codes, and confirmed that the results were in good agreement 
with each other. One of them, described in Gonzalez-Alfonso & Cernicharo (1997, 1999), 
has been recently used to model the far-infrared spectrum of Arp 220 (Gonzalez-Alfonso et 
al. 2004). The other is a Monte Carlo radiative transfer code used as part of a detailed 
study of all the OH lines observed by ISO in galaxies (Smith 2004; Smith et al. 2004). The 
code was developed for the Submillimeter Wave Astronomy Satellite (SWAS) mission and is 
a modified and extended version of the Bernes code (Bernes 1979) but which includes dust 
as well as gas in the radiative transfer, and also corrects some optical depth calculations 
from the original code (Ashby et al. 2000). Both methods are non-local, non-LTE, assume 
spherical symmetry, and include a treatment of continuum photons from dust mixed in with 
the gas. Also, both codes take input as a series of concentric shells, each of which is assigned 
a size, gas and dust temperature, H 2 density, velocity and turbulent velocity width, and 
molecular abundance relative to H 2 . The statistical equilibrium populations of OH in each 
spherical shell are computed by including the excitation by dust emission, excitation through 
collisions, and effects of line trapping. We ran two models to simulate the nucleus of the 
galaxy and the starburst extended ring, described in §5.4.1 and 5.4.2. 
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5. 4-1- Models for the nuclear emission: constraints on the spatial origin of the 79 and 163 

pm OH lines 

We present models for the nuclear OH emission that implicitly assume that the 119 pm 
emission line arises from the nuclear region: X(OH)= 10 -5 is adopted, as well as densities 
~ 5 x 10 5 cm -3 for the bulk of the emitting gas. By assuming a pure nuclear origin for the 
119 pm line, we check whether the other two OH lines could, in such a case, arise from the 
same nuclear region or require a more extended spatial origin. 

The models use the dust parameters derived for the nuclear molecular region described 
in §4. The possible contribution to the OH excitation of far-infrared photons arising from 
the components 1 & 2 of model AGN A is ignored, i.e. only the dust coexistent with the 
molecular gas is taken into account, with predicted flux of 24 — 27 Jy at 35-53 pm (Fig. 4). 

The gas temperature is assumed uniform and equal to 70 K (e.g. Tacconi et al. 1994). 

The H 2 densities derived from the dust model (i.e. a peak density of < n/(H 2 ) > = 

500 cm -3 at the inner radius) are not compatible with the densities inferred from different 
molecular tracers. This indicates that the medium is extremely clumped, as has been also 
argued elsewhere (e.g. Cameron et al. 1993; Tacconi et al. 1994). In order to account ap- 
proximately for this dumpiness in our models, the following strategy is adopted: we use the 
“real” n(H 2 ) ~ 5 x 10 5 cm -3 values for the bulk of the gas, and compute the volume filling 
factor f v =< n(H 2 ) > /n(H 2 ), where the average value is that inferred from the dust model. 

The expected abundances of OH and the dust relative to H 2 , AT (OH) and X(dust), are then 
multiplied by f v , so that the right OH and dust column densities are used in the calculations 
together with the right density values. The same density profile r -1 that was used in the 
dust model is adopted, so that /„ is uniform throughout the nuclear region. 

The modeled fluxes are convolved with the ISO-LWS grating resolution and are com- 
pared with the data in Fig. 6. Solid black lines show the results for the nuclear model that 
assumes X(OH)=10~ 5 and f v = 2 x 10 -5 , the latter value implying a density in the outer 
regions (where the bulk of the emission is generated) of 5 x 10 5 cm -3 . The radial OH column 
density is N(OR) = 2 x 10 17 cm -2 . Besides the 119 pm line, the model nearly reproduces 
the emission in the 79 and 163 pm lines and is consistent with the upper limits given in 
Table 5. 

We also checked the excitation mechanism of the other OH lines by generating an 
additional model with the same parameters as above except for the continuum emission, 
which is now turned off. In this model, therefore, the lines are excited exclusively through 
collisions with H 2 . The resulting flux of the 119 pm line remains unchanged when the dust 
emission is ignored, confirming that the line is collisionally excited. On the other hand, the 
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flux densities of the 79 and 163 /im lines decrease in the “pure-collisional” model by factors 
of 2 and 6, respectively, showing that the emission in these lines is much more affected by 
radiative pumping. We conclude that, if the OH abundance in the nucleus were high enough 
to account for the collisionally excited 119 gm line, the observed fluxes in the 79 and 163 
gm lines can also be explained as arising in the same nuclear region. 


5. 4-2. Models for the starburst emission 

Two simple different approaches have been used to model the OH emission from the 
starburst. First, we have roughly modeled the whole starburst region as a spherical shell 
with external radius of 1.5 kpc, thickness of 0.3 kpc, and average H 2 density < n/(H 2 ) > 
= 3.7 cm -3 , so that the corresponding continuum emission is reproduced with Td = 34 
K (section 4). As shown in section 5.3, the OH 119 /im emission requires densities of 
a few x 10 5 cm -3 , so that we have assumed a volume filling factor f v = 7.5 x 10 -6 and 
therefore a “real” density n(H 2 ) = 5 x 10 5 cm -3 . The kinetic temperature is assumed to 
be Tk = 100 K, and the OH abundance is A(OH)= 2 x 10 -6 x //,, where fh = 0.05 is the 
assumed fraction of warm gas. The result of this model is shown in Fig. 6 (upper dotted 
lines). The 119 /im line is reproduced, but the flux densities of both the 79 gm and 163 /mi 
lines are strongly underestimated. The reason is that the model implicitly assumes that the 
continuum emission, responsible for the excitation of those lines, arises from a very large 
volume, so that the radiation density is weak and has negligible effect on the line excitation. 

Since the OH emission is expected to arise from compact, discrete PDRs in the vicinity 
of O or early B stars, where the continuum infrared radiation density is expected to be 
stronger than assumed above, we have also tried an alternative approach, which consists 
of modelling an individual “typical” cloud of the starburst. We first model the continuum 
from an individual cloud by assuming a central heating source and computing the dust 
equilibrium temperatures at each radial position that result from the balance of heating 
and cooling. The continuum model is adopted if ( i ) leaving aside a scaling factor (N c , the 
number of clouds in the ensemble), the resulting SED is similar to that of the far- infrared 
emission of the starburst (i.e. the 34 K component found in section 4); the value of N c is 
determined by requiring that the absolute continuum flux from the ensemble of clouds is 
equal to that observed for the 34 K component; and ( ii ) we require that the total mass of 
the ensemble does not exceed the mass inferred from the non-nuclear region (< 4 x 10 9 M@). 
Once the continuum is fitted, calculations for OH are performed by assuming Tk =Td, and 
X(OH)= 2 x 10" 6 . 

Several models with various density profiles were found to match the above two require- 
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ments. The common characteristic of all of them is the relatively high column density of 
the individual clouds, N(H. 2 ) > 10 23 cm -2 , which is a consequence of the low effective dust 
temperature (34 K) of the SED. The results of the simplest model, characterized by a flat 
density profile, n(H 2 ) = 5 x 10 5 cm -3 , are given here for reference. With a radius of 5 x 10 17 
cm, a stellar luminosity of 2 x 10 4 L 0 , and N c = 4 x 10 6 , the resulting SED is similar to 
that of the 34 K component (§4). For these clouds we obtain a total mass of 1.8 x 10 9 M©. 
The predicted OH emission/absorption is shown in Fig. 6 (lower dashed lines). In spite of 
the relatively high density and temperatures (28-250 K) throughout the cloud, the 119 pm 
line is predicted to be too weak, and the 79 pm line is predicted in absorption. We have 
found this result quite general: in models where the OH abundance is high enough and the 
radiation density becomes strong enough to pump the 163 pm emission, the continuum at 
79 and 119 pm is absorbed by OH and the predicted emission in the corresponding lines 
is reduced. Models that assume a density profile of r -1 generally predict the 119 pm line 
in absorption. In some models where the OH abundance was allowed to vary with radial 
position, the 79 pm line was predicted in emission but by far too weak to account for the 
observed flux density. 

In conclusion, no starburst model is found to reproduce satisfactorily the emission ob- 
served in the three OH lines. If the local infrared radiation density is strong enough to 
pump the 163 pm line, the other two OH lines are expected to be weak or in absorption. 
Furthermore, the high density assumed for the starburst region would produce a relatively 
high HCN/CO intensity ratio, which is on the contrary ~ 0.01 in the spiral arms (Heifer & 
Blitz 1995). Finally, the PDR models described in §3.2 indicate a density of 1 — 1.5 x 10 3 
cm -3 , i.e. a density much lower than that required to account for the flux density of the 
OH 119 pm line. Therefore, and despite the simplicity of our models, taken together the 
analysis of the OH lines and the derived PDR parameters indicate that the bulk of the OH 
emission arises from X-ray dominated nuclear regions. 


6. CONCLUSIONS 

The main results of this article can be summarized as follows: 

• The complete far-infrared (50-200pm) spectrum of NGC 1068 has been observed for 
the first time. The far-infrared ISO-LWS spectrum has been complemented with the 
mid-infrared data of ISO-SWS and with shorter wavelength (UV, optical and near- 
IR) data from the literature to assemble a composite atomic spectrum as complete as 
possible with the aim of modeling the different line emission components at work. This 
approach has been necessary especially because of the poor spatial resolution of the 
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ISO spectrometers, which were not able to spatially separate the emission components. 

The lines have been interpreted as arising from two physically distinct components: 
the AGN component and a starburst component, the first one nuclear and the second 
one located in the ring at a radius of 15-16 " from the nucleus. Both components are 
characterized by the presence of dust grains, producing strong continuum emission in 
the mid- and far-infrared. The density and ionization parameter of the ~ 5 x 10 6 year 
old starburst are found to be n h ~ 100 cm -3 and log U = -3.5, respectively. Three 
composite models have been computed with different AGN components: the first one 
has the ionizing continuum as derived from Alexander et al. (2000), showing a deep 
trough at energies of a few Rydberg; the second has the monotonically decreasing 
ionizing continuum given by Pier et al. (1994) and the third has a “big blue bump”. 
Two values of the visual extinction (Eg-v = 0.2 and 0.4) have been adopted to correct 
the optical and ultraviolet line fluxes for the reddening. The agreement between the 
composite model with an AGN ionizing continuum characterized by the deep trough 
suggested by Alexander et al. (2000) is very satisfactory, taking into account both 
the simplicity of the photoionization models chosen to avoid dealing with too many 
free parameters and the large number of lines which originate in different physical 
regimes. The agreement between the observed spectrum and what is predicted using 
the canonical ionizing continuum is slightly poorer, while the presence of a big blue 
bump is ruled out. 

• The 50-200pm continuum has been modeled using different components arising from 
both the nucleus and the starburst ring. For the nucleus, we have combined the dust 
emission from the ionized components in the narrow line regions modeled by CLOUDY 
with the neutral component reproduced by the radiative transfer code used for the OH 
molecular emission. For the starburst ring, our CLOUDY modelling of the ionized + 
PDR components could not reproduce the far-infrared emission, while similtaneously 
fitting the far-IR lines. Instead we fit the observed continuum by a neutral molecular 
component, reproduced by two gray body components at temperatures of 20K and 
34K, assuming a steep (/3=2) dust emissivity law. 

• The unique OH emission in the 119 pm line cannot be explained in terms of OH 
excitation through absorption of 35 and 53 pm photons emitted by dust, but rather 
it is collisionally excited. This indicates the presence of a warm and dense region 
with high OH abundance. A simple excitation analysis yields two main alternatives 
for the spatial origin of the observed 119 pm line emission: (i) the nuclear region, 
with 2 x 10 7 M© of warm gas (80 K), an average density of n(H 2 ) = 5 x 10 5 cm -3 , 
and an OH abundance of ~ 10 -5 ; (ii) the starburst region, if ~ 5% of the associated 
mass (~ 6 x 10 7 M©) is warm (~ 100 K), dense (a few x 10 5 cm -3 ), and rich in OH 
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(Y(OH)~ 2 x 10- 6 ). 

• Radiative transfer models that simulate the emission/absorption in all the OH lines 
have been performed for both the nuclear and the starburst region. The models for 
the nucleus quantitatively account for the emission in the three OH lines if the nuclear 
physical conditions pointed out above are assumed. On the other hand, no starburst 
model is found to match the three OH lines simultaneously, because the strong far- 
infrared continuum tends to produce absorption, or to weaken the emission, in the OH 
119 and 79 pm lines (as observed in other galaxies). Therefore, although some con- 
tribution from the extended starburst cannot be ruled out, our models indicate that 
the bulk of the OH emission arises in the nuclear region. The high nuclear OH abun- 
dance required to explain the emission strongly suggest a chemistry deeply influenced 
by X-rays, i.e., an X-ray dominated region. 
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E(Ryd) 


Fig. 1. — The AGN ionizing continua used as input for the photoionization models of 
NGC 1068. The three continua differ in the frequency region between l>E.Rj,d>100, while 
outside this region the Pier et al. (1994) spectrum was adopted. The solid line shows the 
continuum derived from Alexander et al. (2000); the dashed line shows a simple power law 
interpolation; the dotted line shows the presence of the predicted ’’big blue bump”. 
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Fig. 2. — The starburst spectral energy distributions used as input for the photoionization 
models of NGC 1068. The two continua are taken from Leitherer et al. (1999) and represent 
a continuous starburst model (solid line) and an instantaneous model (broken line), both 
with ages of 5 Myr. 
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10 100 
ionization potential I(Ev) 


Fig. 3. — The comparison of the composite models with the observations is shown as the 
ratio of modeled to observed flux ratio for each line, with the ionization potential in the 
x-axis. The assumed reddening is E(B-V)=0.2. Panels from top to bottom: model CM1, 
CM2 and CM3. The short dashed lines represent flux ratios within a factor 3 either ways. 
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Fig. 4. — a) Spectral energy distribution of NGC 1068 and model fit. ISO-SWS fluxes are 
taken from Lutz et al. (2000). The model fit is the composition of ( i ) the NLR components 
1 & 2 of model AGN A, (ii) the emission from the molecular nuclear region, (in) the 34 K 
starburst emission, and (iv) the cold 20 K component, b) and c) Dust temperature versus 
the radial position and radial continuum opacity versus wavelength for the nuclear molecular 
region. 
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Fig. 5. — Energy level diagram of OH. Rotational levels with energies up to 600 K are shown; 
the three lines detected in NGC 1068 are indicated with arrows, as well as the 35 and 53 
/im lines that could play an important role in the radiative excitation. The wavelengths are 
indicated in yum. A-doubling is ignored because the A-doublets are not resolved with the 
ISO grating resolution. 
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Fig. 6. — Comparison between the observed OH lines and model results. As indicated in 
the upper panel, the upper modeled spectrum (solid lines) corresponds to the model for the 
nucleus, the middle one (dotted lines) corresponds to the starburst modelled as a whole, 
and the lower one (dashed lines) corresponds to the starburst modelled as an ensemble of 
individual clouds (see text for details). 
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Table 1. Measured line fluxes from the LWS and SWS grating spectra, with la 

uncertainties. 


Line 

A 

Flux 

Aperture 

reference 


(pm) 

(10~ 13 erg s~ 1 cm~ 2 ) 

r 2 ) 


[Si IX] 3 P 2 - 3 Pi 

2.584 

3.0 

14 x 20 

1 

[Mg VIII] 2 P 3/2 - 2 P 1/2 

3.028 

11. ± 1.1 

14 x 20 

1 

[Si IX] 3 Pi - 3 P 0 

3.936 

5.0 ± 0.6 

14 x 20 

1 

[Mg IV] 2 P 1/2 - 2 P 3/2 

4.487 

7.6 ± 1.5 

14 x 20 

1 

[Ar VI] 2 P 3/2 - 2 P 1/2 

4.529 

15. ± 3. 

14 x 20 

1 

[Fe II] “ 4 F 9/2 - a6 A>/2 

5.340 

5.0 

14 x 20 

1 

[Mg VII] 3 P 2 — ♦ 3 Pi 

5.503 

13. 

14 x 20 

1 

[Mg V] 3 Pi - 3 P 2 

5.610 

18. ± 2. 

14 x 20 

1 

[Ar II] 2 P i/ 2 - 2 P 3/2 

6.985 

13. 

14 x 20 

1 

[Na III] 2 P i /2 - 2 P 3/2 

7.318 

5.8 

14 x 20 

1 

[Ne VI] 2 P 3/2 - 2 P 1/2 

7.652 

110. ± 11. 

14 x 20 

1 

[Fe VII] 3 F, - 3 F 3 

7.815 

3.0 

14 x 20 

1 

[Ar V] 3 P 2 — 3 Pi 

7.902 

< 12. 

14 x 20 

1 

[Na VI] 3 P 2 — 3 Pi 

8.611 

< 16. 

14 x 20 

1 

[Ar III] 3 Pi -> 3 P 2 

8.991 

23.0 ± 3.3 

14 x 20 

1 

[Fe VII] 3 F 3 -» 3 F 2 

9.527 

4.0 

14 x 20 

1 

[S IV] 2 p 3/2 2 p 1/2 

10.510 

58. ± 6. 

14 x 20 

1 

[Ne II] 2 P 3/2 - 2 P 1/2 

12.813 

70. 

14 x 27 

1 

[Ar V] 3 Pi - 3 P 0 

13.102 

< 16. 

14 x 27 

1 

[Ne V] 3 P 2 — 3 P, 

14.322 

97. ± 9.7 

14 x 27 

1 

[Ne III] 3 Pi — 3 P 2 

15.555 

160. ± 32. 

14 x 27 

1 

[Fe II] ° 4 F 7/2 - “ 4 F 9/2 

17.936 

< 10. 

14 x 27 

1 

[S III] 3 P 2 — 3 Pi 

18.713 

40. 

14 x 27 

1 

[Ne V] 3 Pi — 3 P 0 

24.317 

70. ± 7. 

14 x 27 

1 

[O IV] 2 P 3/2 - 2 P 1/2 

25.890 

190. ± 20. 

14 x 27 

1 

[Fe II] “ 6 £> 7/2 - “ 6 D 9/2 

25.988 

8. 

14 x 27 

1 

[S III] 3 Pj -* 3 P 0 

33.481 

55. 

20 x 33 

1 

[Si II] 2 P 3/2 - 2 P l/2 

34.814 

91. 

20 x 33 

1 

[Ne III] 3 P 0 — 3 Pi 

36.013 

18. 

20 x 33 

1 

[O III] 3 P 2 — 3 Pi 

51.814 

114. ± 3. 

80 

2 

[N III] 2 P 3/2 - 2 P 1/2 

57.317 

51.4 ± 2.5 

80 

2 

[O I] 3 Pi - 3 P 2 

63.184 

156. ± 1. 

80 

2 

[O III] 3 P, — 3 P 0 

88.356 

111. ± 1. 

80 

2 

[N II] 3 P 2 — 3 Pi 

121.897 

30.5 ± 1.1 

80 

2 

[O I] 3 P 0 —♦ 3 Pi 

145.525 

11.9 ± 0.4 

80 

2 

[C II] 2 P 3/2 - 2 P 1/2 

157.741 

216. ± 1. 

80 

2 

OH 2 n 1/2 5/2- 2 n 3/2 3/2 

34.60/34.63 

< 3. 

20 x 33 

2 

oh 2 n 1/2 i/2- 2 n 1/2 3/2 

79.11/79.18 

14.4 ± 1.5 

80 

2 

OH 2 n 3/2 5/2- 2 n 3/2 3/2 

119.23/119.44 

11.9 ± 1.2 

80 

2 

oh 2 n 1/2 3/2- 2 n 1/2 i/2 

163.12/163.40 

7.42 ± 0.65 

80 

2 


Note. — (1): from Lutz et al. (2000) and, where errors are available, Alexander et al. (2000); 
(2): this work 
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Table 2. Comparison of observed line fluxes with AGN model predictions 


Line id.A 
(pm ) 

Observed/D 1 /D 2 

Flux (10- 
AGN A model 3 
Comp. 1 Comp. 2 

13 erg s" 1 cm 
AGN B 
Comp. 1 

- 2 ) 

model 4 
Comp. 2 

AGN C model 5 
Comp. 1 Comp. 2 

O VI A .1032+. 1037 

37. 4/4334. /402. 

32.6+19.8 

2.28+1.94 

8.26+5.3 6 

0.31+0.27 

56.7+31.2 

3.41+2.10 

(Lyor)n A .1215 

101.8/3562. /602 

179 

619 

149. 

500 

239. 

638. 

N IV] A .1487 

5.1/103./22.9 

25.2 

6.64 

36.3 

6.28 

59.6 

47.5 

(CIV) n A .1549 

39.7/790. /177. 

142. 

77.6 

110.9 

56.5 

286. 

329. 

Hell A .1640 

21.4/426. /95. 5 

112. 

85.9 

108 

81.8 

346. 

345. 

[Ne V] A .3426 

15.7/95./38.7 

97.7 

23.3 

68.2 

7.36 

271. 

115 

(Ne III] A .3869+. 3968 

19.2/97./43.2 

37.1+11.2 

41.5+12.5 

102.8+31.0 

102. +30.7 

81.6+24.6 

143.+43.0 

Hell A .4686 

6.1/27.6/13. 

15.0 

12.3 

14.7 

11.8 

46.5 

48.4 

(O III] A .4959+. 5007 

2S6./964./496 

86.9+262 

99.5+294 

197.+593 

185.+557. 

148.+446. 

241+726 

[Si VI] A 1.96 

8.0/9.2/8.6 

11.2 

1.53 

10.8 

1.86 

23.3 

3.74 

[Si VII] A 2.48 

8.3 

6.41 

0.1 

9.12 

0.18 

11.0 

0.22 

[Si IX] A 2.584 

3.0 

0.49 

— 

0.84 

— 

0.32 

— 

|Mg VIII] A 3.028 

11. 

2.97 

— 

4.80 

— 

4.29 

— 

(Si IX] A 3.936 

5.4 

0.89 

— 

1.56 

— 

0.56 

— 

[Mg IV] A 4.487 

7.6 

3.85 

7.98 

4.29 

6.92 

9.80 

19.5 

[Ar VI] A 4.529 

15. 

10.3 

4.0 

13.4 

3.10 

26.2 

20.1 

(Mg VII] A 5.503 

13. 

9.68 

0.13 

11.4 

.075 

20.0 

0.45 

[Mg V] A 5.610 

18. 

10.7 

8.63 

7.85 

3.84 

34.6 

29.5 

[Ar II] A 6.985 

13. 

4.56 

7.53 

0.93 

1.36 

0.87 

1.63 

[Na III] A 7.318 

5.8 

0.56 

0.59 

1.27 

1.51 

0.67 

1.21 

[Ne VI] A 7.652 

no. 

153.7 

11.9 

105.4 

3.2 

393. 

62.4 

[Fe VII] A 7.815 

3.0 

1.83 

1.63 

1.65 

0.71 

6.34 

6.30 

[Ar V) A 7.902 

< 12. 

2.14 

1.83 

2.62 

1.92 

3.72 

6.52 

[Na VI] A 8.611 

< 16. 

1.20 

0.16 

0.95 

— 

2.95 

0.65 

[Ar III] + [Mg VII] A 8.991 

25. 

4.63+11.7 

7.34+0.17 

5.97+14.1 

12.4+0.10 

4.83+23.8 

10.7+0.58 

[Fe VII] A 9.527 

4.0 

7.50 

7.0 

6.75 

3.08 

25.8 

26.6 

(S IV] A 10.510 

58. 

38.4 

60.5 

85.3 

126. 

75.8 

187. 

(Ne II] A 12.813 

70. 

5.86 

20.8 

2.20 

4.44 

1.62 

3.26 

[Ar V] A 13.102 

< 16. 

2.23 

2.73 

2.75 

2.88 

3.93 

9.57 

[Ne V] A 14.322 

97. 

91.4 

66.8 

83.2 

28.7 

270. 

273. 

[Ne III] A 15.555 

160. 

44.4 

52.0 

76. 

no. 

42.0 

86.9 

[S III] A 18.713 

40. 

25.8 

56.0 

21.7 

75.7 

18.4 

64.5 

[Ne V] A 24.317 

70. 

40.8 

52.6 

35.8 

22.1 

116. 

215. 

(O IV] A 25.890 

190. 

24.1 

71.9 

28.3 

69.7 

56.4 

240. 

(S III] A 33.481 

55. 

8.00 

29.8 

5.57 

41.1 

4.84 

38.4 

[Si II] A 34.814 

91. 

12.9 

22.6 

6.01 

24.2 

5.28 

29.4 

[Ne III] A 36.013 

18. 

3.32 

4.4 

5.74 

9.40 

3.20 

7.46 

[O III] A 51.814 

110. 

9.04 

33.9 

16.3 

63.2 

9.06 

48.5 

[N III] A 57.317 

51. 

2.23 

11.4 

2.83 

18.9 

1.88 

16.1 

(O I] A 63.184 

156. 

3.84 

1.8 

1.57 

2.26 

1.24 

3.10 

[O III] A 88.356 

110. 

1.22 

8.2 

2.25 

15.5 

1.23 

12.6 

[N II] A 121.897 

30. 

0.28 

0.87 

.07 

0.40 

.06 

0.52 

[O I] A 145.525 

12. 

0.24 

0.14 

.09 

0.18 

.07 

0.24 

[C II] A 157.741 

220. 

0.50 

1.58 

0.18 

0.97 

.14 

1.09 


This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article 
may differ from this preprint. 

Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/ , 428495 M . 

- 33 - 


Table 2 — Continued 


Line id. A Flux (10 13 erg s 1 cm 2 ) 

(pm ) Observed/D 1 /D 2 AGN A model 3 AGN B model 4 AGN C model 5 

Comp. 1 Comp. 2 Comp. 1 Comp. 2 Comp. 1 Comp. 2 


1 Dereddened line flux, assuming E b-V = 0*4 

2 Dereddened line flux, assuming Eg-v =0.2 

3 AGN A parameters: component 1: Log U=-l., Log n=4, internal radius ~ 21 pc, external radius ~ 
109 pc, ionizing spectrum from Alexander et al. (2000): component 2: Log U=-2., Log n=3.3, internal 
radius ~ 153 pc, external radius ~ 362 pc, ionizing spectrum from Alexander et al. (2000). 

4 AGN B parameters: same as AGN A models, but with the ionizing spectrum from Pier et al. (1994) 

5 AGN C parameters: same as AGN A models, but with the ionizing spectrum that includes a big blue 
bump (see text) 
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Table 3. Comparison of observed line fluxes with the Ring Starburst model predictions 


Line id. A Flux (10“ 13 erg s” 1 cm"" 2 ) 


(pm ) 

Observed 

SBR A 1 

SBR B 2 

SBR C 3 

SBR D 4 

SBR E 5 

SBR F 6 

O VI A .1032+. 1037 

37.4/4334. /402. 

— 

— 

— 







(Lya)n A .1215 

101. 8/3562. /602 

1330. 

865. 

997. 

875. 

1010. 

917. 

N IV] A .1487 

5.1/103. /22. 9 

71.9 

0.20 

85.8 

.020 

93.7 

0.03 

(CIV)n A .1549 

39. 7/790./ 177. 

795. 

0.70 

858. 

0.80 

924. 

0.87 

Hell A .1640 

21.4/426./95.5 

693. 

50.2 

696. 

50.5 

696. 

50.8 

[Ne V] A .3426 

15.7/95. /38. 7 

58.1 

— 

69.6 

— 

73.9 

— 

[Ne III] A .3869+. 3968 

19.2/97./43.2 

314.+94.7 

31.1 +9.4 

309.493.1 

31.549.5 

316.495.7 

32.49.7 

Hell A .4686 

6.1/27.6/13. 

97.3 

7.4 

97.7 

7.4 

97.3 

7.4 

[O III] A .4959+. 5007 

2S6./964./496 

1330.4-3990. 

25.1475.6 

138344158. 

26.0478.2 

1412.44257. 

26.8480.5 

[Si VI] A 1.96 

8. 0/9. 2/8. 6 

— 

— 

— 

— 

— 

— 

[Si VII] A 2.48 

8.3 

— 

— 

— 

— 

— 

— 

[Si IX] A 2.584 

3.0 

— 

— 

— 

— 

— 

— 

[Mg VIII] A 3.028 

11. 

— 

— 

— 

— 

— 

— 

[Si IX] A 3.936 

5.4 

— 

— 

— 

— 

— 

— 

[Mg IV] A 4.487 

7.6 

3.30 

— 

3.50 

— 

3.53 

— 

[Ar VI] A 4.529 

15. 

6.20 

— 

7.72 

— 

8.02 

— 

[Mg VII] A 5.503 

13. 

— 

__ 

— 

— 

— 

— 

[Mg V] A 5.610 

18. 

3.10 

— 

3.63 

— 

3.73 

— 

[Ar II] A 6.985 

13. 

5.90 

3.70 

4.29 

3.70 

4.12 

3.73 

[Na III] A 7.318 

5.8 

0.60 

0.10 

0.54 

0.11 

0.54 

0.11 

[Ne VI] A 7.652 

110. 

16.9 

— 

21.1 

— 

21.8 

— 

[Fe VII] A 7.815 

3.0 

0.40 

— 

0.48 

— 

0.50 

— 

[Ar V] A 7.902 

<12. 

7.50 

— 

8.84 

— 

9.11 

— 

[Na VI] A 8.611 

<16. 

— 

— 

— 

— 

— 

— 

[Ar III] + [Mg VII] A 8.991 

25. 

46.5 

7.10 

43.2 

7.13 

42.2 

7.13 

[Fe VII] A 9.527 

4.0 

1.70 

— 

2.04 

— 

2.11 

— 

[S IV] A 10.510 

58. 

312. 

0.80 

346. 

0.84 

333. 

0.81 

[Ne II] A 12.813 

70. 

7.80 

5.50 

5.44 

5.38 

5.38 

5.61 

[Ar V] A 13.102 

<16. 

21.9 

— 

25.5 

— 

24.9 

— 

[Ne V] A 14.322 

97. 

271. 

— 

315. 

— 

314. 

— 

[Ne III] A 15.555 

160. 

295. 

50.5 

276. 

50.2 

270. 

49.8 

[S III] A 18.713 

40. 

180. 

19.2 

178. 

20.3 

203. 

23.7 

[Ne V] A 24.317 

70. 

303. 

— 

346. 

— 

294. 

— 

[OIV] A 25.890 

190. 

1660. 

1.90 

1762. 

2.00 

1419. 

1.62 

[S III] A 33.481 

55. 

785. 

87.5 

673. 

81.2 

350. 

46.2 

[Si II] A 34.814 

91. 

126. 

54.8 

89.4 

49.5 

41.9 

29.4 

[Ne III] A 36.013 

18. 

52.8 

8.90 

49.2 

8.91 

47.5 

8.74 

[O III] A 51.814 

110. 

2980. 

81.5 

3217 

87.5 

2524. 

71.9 

[N III] A 57.317 

51. 

1110. 

40.0 

921 

35.0 

372.9 

14.7 

[O I] A 63.184 

156. 

165. 

130. 

113. 

127. 

116. 

124. 

[O III] A 88.356 

110. 

4980. 

138. 

3729. 

109. 

977. 

29.0 

[N II] A 121.897 

30. 

80.0 

37.3 

47.8 

29.7 

12.3 

8.78 

[O I] A 145.525 

12. 

16.2 

12.9 

11.0 

12.5 

10.1 

11.3 

[C II] A 157.741 

220. 

739. 

465. 

181. 

172. 

67.6 

46.5 



This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article 
may differ from this preprint. 

Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi: 10. 1086/' 428495". 

- 35 - 


Note. — *: this line was used for normalization 

1 SBR A parameters: Log U=-2.5, Log n=1.0, ionizing spectrum from Starburst99 with instantaneous star-formation law, M = 
10 6 Af©, IMF: =2.35 Mttp = 100 Af©, Mf ou , = 1A/©, nebular emission included, Z=0.020, age of 5 Myr. The integration was stopped at 
a temperature of 50K, the adopted abundances are those relative to HII regions and grain emission is included. The adopted number 
of clouds is 33000. 

2 SBR B parameters: Log U=-3.5, Log n=1.0, all other parameters as for SBR A. The number of clouds adopted is 33000 
3 SBR C parameters: Log U=-2.5, Log n=2.0, all other parameters as for SBR A. The number of clouds adopted is 3300 

4 SBR D parameters: Log U=-3.5, Log n=2.0, all other parameters as for SBR A. The number of clouds adopted is 3300 

5 SBR E parameters: Log U=-2.5, Log n=3.0, all other parameters as for SBR A. The number of clouds adopted is 330 

6 SBR F parameters: Log U=-3.5, Log n=3.0, all other parameters as for SBR A. The number of clouds adopted is 330 
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Table 4. Comparison of observed line fluxes with composite model predictions 


Line id.A 

(A* m ) 

Flux (10- 13 
Observed/D 1 /D 2 

erg s 1 
CM1 3 

cm 2 ) 
CM2 4 

CM3 5 

0 VI A .1032+. 1037 

37.4/4334./402. 

56.6 

14.2 

93.4 

(Lya)n A .1215 

101.8/3562./602 

1673. 

1524. 

1752. 

N IV] A .1487 

5.1/103./22.9 

31.8 

42.6 

107. 

(CIV)n A .1549 

39.7/790./177. 

220. 

168. 

616. 

Hell A .1640 

21.4/426./95.5 

248. 

240. 

741. 

[Ne V] A .3426 

15.7/95./38.7 

121. 

75.6 

386. 

[Ne III] A .3869+.3968 

19.2/97./43.2 

143. 

307. 

333. 

Hell A .4686 

6.1/27.6/13. 

34.7 

33.9 

102.3 

[O III] A .4959+. 5007 

256./964./496 

847. 

1636. 

1665. 

[Si VI] A 1.96 

8.0/9.2/8.6 

12.7 

12.7 

27.0 

[Si VII] A 2.48 

8.3 

6.5 

9.3 

11.2 

[Si IX] A 2.584 

3.0 

0.5 

0.8 

0.3 

[Mg VIII] A 3.028 

11. 

3.0 

4.8 

4.3 

[Si IX] A 3.936 

5.4 

0.9 

1.6 

0.6 

[Mg IV] A 4.487 

7.6 

11.8 

11.2 

29.3 

[Ar VI] A 4.529 

15. 

14.3 

16.5 

46.3 

[Mg VII] A 5.503 

13. 

9.8 

11.5 

20.5 

[Mg V] A 5.610 

18. 

19.3 

11.7 

64.1 

[Ar II] A 6.985 

13. 

15.8 

6.0 

6.2 

[Na III] A 7.318 

5.8 

1.3 

2.8 

2.0 

(Ne VI] A 7.652 

110. 

166. 

109. 

455. 

[Fe VII] A 7.815 

3.0 

3.5 

2.4 

12.6 

[Ar V] A 7.902 

<12. 

4.0 

4.5 

10.2 

[Na VI] A 8.611 

<16. 

1.4 

0.9 

3.6 

[Ar III] + [Mg VII] A 8.991 

25. 

31. 

39.7 

47.0 

[Fe VII] A 9.527 

4.0 

14.5 

9.8 

52.4 

[S IV] A 10.510 

58. 

99.7 

212. 

263. 

[Ne II) A 12.813 

70. 

32.0 

12.0 

10.3 

[Ar V] A 13.102 

<16. 

5.0 

5.6 

13.5 

[Ne V] A 14.322 

97. 

158. 

112. 

543. 

[Ne III] A 15.555 

160. 

147. 

236. 

179. 

[S III] A 18.713 

40. 

102. 

118. 

103. 

[Ne V] A 24.317 

70. 

93.4 

57.9 

331. 

[O IV] A 25.890 

190. 

98. 

100. 

298. 

[S III] A 33.481 

55. 

119. 

128. 

124. 

[Si II] A 34.814 

91. 

85. 

79.7 

84.2 

[Ne III] A 36.013 

18. 

16.6 

24.0 

19.6 

[O III] A 51.814 

no. 

130. 

167. 

145. 

[N III] A 57.317 

51. 

48.6 

56.7 

53. 

[O I] A 63.184 

156. 

133. 

131. 

131. 

[O III] A 88.356 

110. 

118. 

127. 

123. 

[N II] A 121.897 

30. 

30.8 

30. 

30.3 

[O I] A 145.525 

12. 

12.9 

12.8 

12.8 

[C II] A 157.741 

220. 

174. 

173. 

173. 


1 Dereddened line flux, assuming E#-v = 0.4 
2 Dereddened line flux, assuming E b-v =0.2 
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Table 5. Comparison of observed line fluxes with model predictions for the nuclear region 


Line id. A 

Flux (10- 12 
Observed 

erg s 1 cm 2 ) 
Modeled 

Notes 

34/xm 

< 0.1 

-0.43 

(absorption) 

53fim 

< 1.2 

-0.66 

(absorption) 

79fxm 

1.1 

1.13 


84/xm 

< 1.2 

0.06 


98/im 

< 1.2 

0.17 


119/im 

1.3 

1.60 


163/xm 

0.38 

0.35 



